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Abstract 

We present an exact one-loop calculation of the tunneling process in Euclidean 
quantum gravity describing creation of black hole pairs in a de Sitter universe. 
Such processes are mediated by S 2 x S 2 gravitational instantons giving an imagi- 
nary contribution to the partition function. The required energy is provided by the 
expansion of the universe. We utilize the thermal properties of de Sitter space to 
describe the process as the decay of a metastable thermal state. Within the Eu- 
clidean path integral approach to gravity, we explicitly determine the spectra of the 
fluctuation operators, exactly calculate the one-loop fluctuation determinants in the 
C-function regularization scheme, and check the agreement with the expected scal- 
ing behaviour. Our results show a constant volume density of created black holes 
at late times, and a very strong suppression of the nucleation rate for small values 
of A. 
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1 Introduction 



Instantons play an important role in flat space gauge field theory [45]. Being 
stationary points of the Euclidean action, they give the dominant contribu- 
tion to the Euclidean path integral thus accounting for a variety of important 
phenomena in QCD-type theories. In addition, self-dual instantons admit su- 
persymmetric extensions, which makes them an important tool for verifying 
various duality conjectures like the AdS/CFT correspondence [4]. More gen- 
erally, the Euclidean approach has become the standard method of quantum 
field theories in flat space. 

Since the theory of gravity and Yang-Mills theory are somewhat similar, it is 
natural to study also gravitational instantons. An impressive amount of work 
has been done in this direction, leading to a number of important discoveries. 
A thorough study of instanton solutions of the vacuum Einstein equations 
and also those with a A-term has been carried out [20,17,24]. These solutions 
dominate the path integral of Euclidean quantum gravity, leading to interest- 
ing phenomena like black hole nucleation and quantum creation of universes. 
Perhaps one of the most spectacular achievements of the Euclidean approach 
is the derivation of black hole entropy from the action of the Schwarzschild in- 
stanton [22] . In addition, gravitational instantons are used in the Kaluza-Klein 
reductions of string theory. 

Along with these very suggestive results, the difficulties of Euclidean quan- 
tum gravity have been revealed. Apart from the usual problem of the non- 
renormalizability of gravity, which can probably be resolved only at the level of 
a more fundamental theory like string theory, the Euclidean approach presents 
other challenging problems. In field theories in flat space the correlation func- 
tions of field operators are holomorphic functions of the global coordinates 
in a domain that includes negative imaginary values of the time coordinate, 
t = —it, where r is real and positive [48]. This allows one to perform the 
analysis in the Euclidean section and then analytically continue the functions 
back to the Lorentzian sector to obtain the physical predictions. In curved 
space the theorems that would ensure the analyticity of any quantities arising 
in quantum gravity are not known. As a result, even if Euclidean calculations 
make sense, it is not in general clear how to relate their result to the Lorentzian 
physics. 

This difficulty is most strikingly illustrated by the famous problem of the con- 
formal sector in Euclidean quantum gravity. If one tries evaluating the path 
integral over Riemannian metrics, then one discovers that it diverges because 
the Euclidean gravitational action is not bounded from below and can be 
made arbitrarily large and negative by conformal rescaling of the metric [25]. 
Such a result is actually expected, for if the integral did converge (with some 
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regularization), then one could give a well-defined meaning to the canonical 
ensemble of the quantum gravitational field. However, the possibility of having 
a black hole causes the canonical ensemble to break down - since the degen- 
eracy of black hole states grows faster than the Boltzmann factor decreases. 
One can, 'improve' the Euclidean gravitational action by analytically contin- 
uing the conformal modes, let us call them h, via h — * ih, and this improves 
the convergence of the integral [25]. This shows that if there is a well-defined 
Euclidean path integral for the gravitational field, then the relation to the 
Lorentzian sector is more complicated than just via t — > —ir. 

Unfortunately, it is unknown at present whether one can in the general case 
find a physically well-defined and convergent path integral for the gravitational 
field. At the same time, the idea of constructing it is conceptually simple [46]: 
one should start from the Hamiltonian path integral over the physical degrees 
of freedom of the gravitational field. Such an integral certainly makes sense 
physically and is well-convergent, since the Hamiltonian is positive - at least 
in the asymptotically flat case. The Hamiltonian approach is not covariant, 
but one can covariantize it by changing the integration variables, which leads 
to a manifestly covariant and convergent path integral for gravity. The main 
problem with this program is that in the general case it is unclear how to isolate 
the physical degrees of freedom of the gravitational field. For this reason, so 
far the program has been carried out only for weak fields in the asymptotically 
flat case [46]. Remarkably, the result has been shown to exactly correspond 
to the the standard Euclidean path integral with the conformal modes being 
complex-rotated via h — > ih. This lends support to the Euclidean approach 
in gravity and allows one to hope that the difficulties of the method can be 
consistently resolved; (see, for example, [7,6] for the recent new developments 
within the lattice approach). 

One can adopt the viewpoint that Euclidean quantum gravity is a meaningful 
theory within its range of applicability, at least at one-loop level, by assum- 
ing that a consistent resolution of its difficulties exists. Then already in its 
present status the theory can be used for calculating certain processes, most 
notably for describing tunneling phenomena, in which case the Euclidean am- 
plitude directly determines the probability. The analytic continuation to the 
Lorentzian sector in this case is not necessary, apart from when the issue of 
the interpretation of the corresponding gravitational instanton is considered. 
The important example of a tunneling process in quantum gravity is the cre- 
ation of black holes in external fields. Black holes are created whenever the 
energy pumped into the system is enough in order to make a pair of virtual 
black holes real [33]. The energy can be provided by the heat bath [30,38,5], by 
the background magnetic field [21,19,16,15], by the expansion of the universe 
[28,10,41], by cosmic strings [37], domain walls [11], etc; (see also [43,35,36]). 
Besides, one can consider pair creation of extended multidimensional objects 
like p-branes due to interaction with the background supergravity fields [14]. 
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In all these examples the process is mediated by the corresponding gravita- 
tional instanton, and the semiclassical nucleation rate for a pair of objects on 
a given background is given by 

r = Aexp{-(/ obj -/ bg )} . (1.1) 

Here J bj is the classical action of the gravitational instanton mediating cre- 
ation of the objects, is the action of the background fields alone, and the 
prefactor A includes quantum corrections. 

In most cases the existing calculations of black hole pair creation processes 
consider only the classical term in (1.1). This is easily understood, since loop 
calculations in quantum gravity for non-trivial backgrounds are extremely 
complicated. To our knowledge, there is only one example of a next-to-leading- 
order computation, which was undertaken in [30] by Gross, Perry, and Yaffe 
for the Schwarzschild instanton background. The aim of the present paper is 
to consider one more example of a complete one-loop computation in quantum 
gravity. 

The problem we are interested in is the quantum creation of black holes in de 
Sitter space. This problem was considered by Ginsparg and Perry [28], who 
identified the instanton responsible for this process, which is the S 2 x S* 2 so- 
lution of the Euclidean Einstein equations R^ u = Kg^ u for A > 0. Ginsparg 
and Perry noticed that this solution has one negative mode in the physi- 
cal sector, which renders the partition function complex, thus indicating the 
quasi-classical instability of the system. This instability leads to spontaneous 
nucleation of black holes in the rapidly inflating universe. This is the domi- 
nant instability of de Sitter space, since classically the space is stable [28] . The 
energy necessary for the nucleation is provided by the A-term, which drives 
different parts of the universe apart thereby drugging the members of a virtual 
black hole pair away from each other. The typical radius of the created black 
holes is 1/y/X, while the the nucleation rate is of the order of exp(— it /KG), 
where G is Newton's constant. As a result, for AG ~ 1 when inflation is fast, 
the black holes are produced in abundance but they are small and presumably 
almost immediately evaporate. Large black holes emerge for AG <C 1 when 
inflation slows down, and these can probably exist for a long time, but the 
probability of their creation is exponentially small. This scenario was further 
studied in Refs. [10,9,18] (see also references in [9]), where the generalization 
to the charged case was considered and also the subsequent evolution of the 
created black holes was analyzed. However, the one-loop contribution so far 
has not been computed. 

A remarkable feature of the S 2 x S 2 instanton is its high symmetry. In what 
follows, we shall utilize this symmetry in order to explicitly determine spectra 
of all relevant fluctuation operators in the problem. We shall use the C-function 
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regularization scheme in order to compute the one-loop determinants, which 
will give us the partition function Z[S 2 x S 2 } for the small fluctuations around 
the S 2 x S 2 instanton. We shall then need to normalize this result. The nor- 
malization coefficients is ^[S* 4 ], the partition function for small fluctuations 
around the S 4 instanton, which is the Euclidean version of the de Sitter space. 
The one-loop quantization around the S 4 instanton was considered by Gib- 
bons and Perry [27], and by Christensen and Duff [13], but unfortunately in 
none of these papers the analysis was completed. We shall therefore reconsider 
the problem by rederiving the spectra of fluctuations around S 4 and comput- 
ing the determinants within the ^-function scheme, thereby obtaining a closed 
one- loop expression for ZfS* 4 ]. 

In our treatment of the path integral we follow the approach of Gibbons and 
Perry [27]; (see also [42]). In order to have control over the results, we work in a 
one-parameter family of covariant gauges and perform the Hodge decomposi- 
tion of the fluctuations. These are then expanded with respect to the complete 
sets of basis harmonics, and the perturbative path integration measure is de- 
fined as the square root of the determinant of the metric on the function space 
of fluctuations. To insure the convergence of the integral over the conformal 
modes, which enter the action with the wrong sign, we essentially follow the 
standard recipe h —> ih [25]; (see also Ref.[42], where a slightly disguised form 
of the same prescription was advocated). The subtle issue is that the confor- 
mal operator Ao = — 3V M V M — 4A has a finite number, A/", of negative modes, 
and these enter the action with the correct sign from the very beginning. 
Our treatment of these special modes is different from that by Hawking [32], 
who suggests that such modes should be complex-rotated twice, the partition 
function then acquiring the overall factor of . However, the presence of this 
factor in the partition function would lead to unsatisfactory results, and on 
these grounds we are led to not rotating the special conformal modes at all. 

The path integral is computed by integrating over the Fourier expansion co- 
efficients, which leads to infinite products over the eigenvalues. The only con- 
formal modes giving contribution to the result are the special negative modes 
discussed above. We carefully analyze the resulting products to make sure that 
all modes are taken into account and that the dependence of the gauge-fixing 
parameter cancels thereby indicating the correctness of the procedure. We 
give a detailed consideration to the zero modes of the Faddeev-Popov opera- 
tor, which arise due to the background isometries. The integration over these 
modes requires a non-perturbative extension of the path-integration measure, 
and we find such a non-perturbative measure in the zero mode sector to be 
proportional to the Haar measure of the isometry group. Collecting all terms 
yields the partition function for small fluctuations around a background in- 
stanton configuration in terms of infinite products over eigenvalues of the 
gauge-invariant operators. We then use the explicitly known spectra of fluc- 
tuations around the S 2 x S* 2 and S 4 backgrounds in order to calculate the 
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partition functions. 



The rest of the paper is organized as follows. In Sec. 2 we present our deriva- 
tion of the black hole nucleation rate within the finite temperature approach. 
In Sec. 3 the path integration procedure is considered. The spectra of small 
fluctuations around the S 2 x S 2 instanton are computed in Sec. 4 via a direct 
solving of the differential equations in the eigenvalue problems. The spectra of 
the fluctuations around the S A instanton are rederived in Sec. 5 with the use 
of group theoretic arguments. The partition functions are computed in Sec. 6, 
and Sec. 7 contains the final expression for the black hole nucleation rate to- 
gether with some remarks. We present a detailed analysis of the ^-functions 
in the Appendix. We use units where c = h — &b = 1- 



2 Black hole nucleation rate 

In this section we shall derive the basic formula for the black hole nucleation 
rate in de Sitter space, whose different parts will be evaluated in the next 
sections. The existing derivations of the nucleation rate [28,10] recover only the 
classical factor in (1.1). In addition, it is not always clear to which volume the 
rate refers. We argue that our formula (2.15) gives the nucleation probability 
per Hubble volume and unit time as measured by a freely falling observer. The 
basic idea of our approach is to utilize the relation between the inflation and 
thermal properties of de Sitter space. This will allow us to use the standard 
theory of decay of metastable thermal states [39,40,3]. 

Let us consider the partition function for the gravitational field 



where the integral is taken over Riemannian metrics, and / = I[g^ v } is the 
Euclidean action for gravity with a positive A terms; see Eq.(3.1) below. The 
path integration procedure will be considered in detail in the next section. At 
present let us only recall that in the semiclassical approximation the integral 
is approximated by the sum over the classical extrema of the action /, that is 



Here Z\ = Z\M.j\ is the partition function for the small gravitational fluctu- 
ations around a background manifold A4i with a metric g l subject to the 




(2.2) 



(2.3) 
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Euclidean Einstein equations = Ag^. Schematically one has 

z exp(-g 



where J/ is the classical action for the Z-th extremum, and A; is the operator 
for the small fluctuations around this background. 

The dominant contribution to the sum in (2.4) is given by the S 4 instanton, 



which is the four-dimensional sphere with the radius y3/A and the standard 
metric. Since this is a maximally symmetry space, its action / = —3tt/AG is 
less than that of any other instanton. Hence, 

Z,Z[5l = ^5. (2.5) 
1 J v/DetA V ; 



On the other hand, the S* 4 instanton describes the thermal properties of de 
Sitter space [22,23], since it can be obtained by an analytic continuation via 
t — > r = it of the region of the de Sitter solution 

ds 2 = -(1 - -r 2 ) dt 2 + -^j— + r 2 {dd 2 + sin 2 dd V 2 ) (2.6) 
3 1 — %r z 



contained inside the event horizon, r < y 3/A. Let us call this region a Hubble 
region. Its boundary, the horizon, has the area A = 127r/A. The temperature 
associated with this horizon is T = the entropy S = A/ 'AG = 3n/AG 

and the free energy F = —TS. The same values can be obtained by writing 
the partition function for the S 14 instanton as 

Z[S A ] = e- pF (2.7) 



with (3 = 1/T. Indeed, since S* 4 is periodic in all four coordinates, any of 
them can be chosen to be the 'imaginary time'. The corresponding period, 
/3 = 27ry^, can be identified with the proper length of a geodesic on S* 4 , 
all of which are circles with the same length. This gives the correct de Sitter 
temperature. Comparing (2.7) and (2.5) one obtains (3F = — 37r/AG+. . . , the 
dots denoting the quantum corrections, and this again agrees with the result 
for the de Sitter space. To recapitulate, the partition function of quantum 
gravity with A > is approximately 

Z w e"^ , (2.8) 
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Fig. 1. The leading contribution to the partition function comes from the S 4 and 
S 2 x S 2 gravitational bubbles, the effect of the latter being purely imaginary. 

where 1//3 is the de Sitter temperature and F is the free energy in the Hubble 
region. 

Let us now consider the contribution of the other instantons. One has 



where the prime indicates that Aii ^ S A . Now, for AG <C 1 all terms in 
the sum are exponentially small and can safely be neglected as compared to 
the unity, if only they are real. If there are complex terms, then they will 
give an exponentially small imaginary contribution. The S 2 x S 2 instanton is 
distinguished by the fact that its partition function is purely imaginary due 
to the negative mode in the physical sector [28]. This is the only solution for 
A > which is not a local minimum of the action in the class of metrics with 
constant scalar curvature [20]. Hence (see Fig.l), 

BF ( Z[S 2 xS 2 ]\ ( / Z[S 2 xS 2 ]\\ , , 



where Z[S 2 x S 2 } is purely imaginary. As a result, the partition function can 
still be represented as Z ~ e _/3F , where the real part of F is the free energy 
of the Hubble region, and the exponentially small imaginary part is given by 

w=-%?- (2 ' n) 



It is natural to relate this imaginary quantity also to the free energy. We are 
therefore led to the conclusion that the free energy of the Hubble region has a 
small imaginary part, thus indicating that the system is metastable. The decay 
of this metastable state will lead to a spontaneous nucleation of a black hole 
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in the Hubble region, which can be inferred from the geometrical properties 
of the S 2 x S 2 instanton. 

The S 2 x S 2 instanton can be obtained via the analytic continuation of the 
Schwarzschild-de Sitter solution [26,28,10] 

dv 2 

ds 2 = -Ndt 2 + — + r 2 {dd 2 + sm 2 dd V 2 ). (2.12) 



Here N = l-^--f r 2 , and for 9M 2 A < 1 this function has roots at r = r + > 
(black hole horizon) and at r = r ++ > r + (cosmological horizon). One has 
N > for r + < r < r ++ , and only this portion of the solution can be 
analytically continued to the Euclidean sector via t — > r = it. The conical 
singularity at either of the horizons can be removed by a suitable identification 
of the imaginary time. However, since the two horizons have different surface 
gravities, the second conical singularity will survive. The situation improves 
in the extreme limit, r + — > r ++ — > -7=, since the surface gravities are then 
the same and both conical singularities can be removed at the same time. 
Although one might think that the Euclidean region shrinks to zero when 
the two horizons merge, this is not so. The limit r + — > r ++ implies that 
9M 2 A = 1 — 3e 2 with e — > 0. One can introduce new coordinates $1 and </?i via 
cos^i = (y/Xr — l)/e + e/6 and (pi = y/Xer. Passing to the new coordinates 
and taking the limit e — > 0, the result is 

ds 2 = \( d^ 2 + sin 2 0i dv?i 2 + dd 2 + sin 2 , (2- 13) 



and this S* 2 x S 2 metric fulfills the Einstein equations. Since the instanton 
field determines the initial value for the created real time configuration, one 
concludes that the S* 2 x S* 2 instanton is responsible for the creation of a black 
hole in the Hubble region. This black hole fills the whole region, since its size 
is equal to the radius of the cosmological horizon. 



It is well known that the region r < ^3/A of the static coordinate system in 
(2.6) covers only a small portion of the de Sitter hyperboloid [47]; (see Fig.2). 
In order to cover the whole space, one can introduce an infinite number of freely 
falling observers and associate the interior of the static coordinate system with 
each of them. Hence, the spacetime contains infinitely many Hubble regions. 
It is also instructive to use global coordinates covering the whole de Sitter 
space, 

ds 2 = (~de + d X 2 + sin 2 X (^ 2 + sin 2 M V 2 )\ , (2.14) 



where £ G [— 7r/2, tt/2] and x £ [0, 7r]. The trajectory of a freely falling observer 
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Fig. 2. Left: The conformal diagram of de Sitter space in coordinates (2.14). 
The trajectories x =cons t- are timelike geodesies. The diamond-shaped region in 
the center is the Hubble region of the geodesic observer at \ = %/2. Although 
this region completely covers the hypersurface So, at later times one needs more 
observers to cover the hypersurface Si with the interiors of their horizons - the 
Hubble regions proliferate. Right: The de Sitter hyperboloid in the embedding 
Minkowski space (with two dimensions suppressed). The Hubble region of the 
inertial observer moving along the hyperbola x = 0, y > is the portion of the 
hyperboloid lying to the right from the two shaded strips. This corresponds to 
the interior region of the observer's static coordinate system. 

is X — Xo ( an d also $ = i? 0) V 9 = <A))> an d the domain of the associated static 
coordinate system, the Hubble region, is the intersection of the interiors of the 
observer's past and future horizons [34]. Let £ be a spacelike hypersurface, 
say £ = £o- If £o = then £ is completely contained inside the Hubble region 
of a single observer with x = 7r/2 (see Fig. 2). However, for late moments of 
time, £ — > 7r/2, one needs more and more independent observers in order to 
completely cover £ by the union of their Hubble regions. One can say that 
the Hubble regions proliferate with the expansion of the universe. 

Since de Sitter space consists of infinitely many Hubble regions, the black hole 
nucleation will lead to some of the regions being completely filled by a black 
hole, but most of the regions will be empty. The number of the filled regions 
divided by the number of those without a black hole is the probability for a 
black hole nucleation in one region. This is proportional to Q(F) in (2.11). 

One can argue that the black holes are actually created in pairs [33,36], 
where the two members of the pair are located at the antipodal points of 
the de Sitter hyperboloid. This can be inferred from the conformal diagram 
of the Schwarzschild-de Sitter solution, which contains an infinite sequence of 
black hole singularities and spacelike infinities; see Fig.3. One can identify the 
asymptotically de Sitter regions in the diagram related by a horizontal shift, 
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and the spacetime will then consist of two black holes at antipodal points 
of the closed universe. This agrees with the standard picture of particles in 
external fields being created in pairs. 



r = r = oo 




r = 



Fig. 3. The conformal digram for the extreme Schwarzschild-de Sitter solution. 

The surface gravity of the extreme Schwarzschild-de Sitter solution is finite 
when defined with respect to the suitably normalized Killing vector [10]. This 
gives a non-zero value for the temperature of the nucleated black holes, which 
can be read off also from the S 2 x S 2 metric: it is the inverse proper length of 
the equator of any of the two spheres, Tbh = % ■ How can it be that this is 
different from the temperature of the heat bath, which is the de Sitter space 
with Tds = ? F° r example, in the hot Minkowski space the nucleated 

black holes have the same temperature as the heat bath [30]. However, the 
global structure of de Sitter space is different from that of Minkowski space. 
The fluctuations cannot absorb energy from and emit energy into the whole of 
de Sitter space, but can only exchange energy with the Hubble region. Thus 
the energy exchange is restricted. As a result, the local temperature in the 
vicinity of a created defect may be different from that of the heat bath, but 
reduces to the latter in the asymptotic region far beyond the cosmological 
horizon. 

The relation of the imaginary part of the free energy to the rate of decay of a 
metastable thermal state Y was considered in [39,40,3]. If the decay is only due 
to tunneling then Y = 23(F). Suppose that there is an additional possibility 
to classically jump over the potential barrier. In this case on top of the barrier 
there is a classical saddle point configuration whose real time decay rate is de- 
termined by the saddle negative mode co>_ . At low temperatures the tunneling 
formula is then still correct, while for T > one has Y = M=r ^s(F). In our 
problem the saddle point configuration also exists, the S 2 x S 2 instanton, but 
its real time analog, the Schwarzschild-de Sitter black hole, is stable. It seems 
therefore that there is no classical contribution to the process and the black 
hole nucleation is a purely quantum phenomenon.^] [One can imagine that the 
effective potential barrier is infinitely high, such that a classical transition is 

2 We do not understand the classical interpretation of the Euclidean saddle point 
solution suggested in [30]. The argument uses a family of non-normalizable de- 
formations of the instanton, and the action is finite as long as they are 'static'. 
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forbidden, but at the same time so narrow that the tunneling rate is non-zero.] 
As a result, the rate of quasiclassical decay of the de Sitter space is given by 
r = 23(F). Using Eq.(2.11), 



Here T = is the temperature of the de Sitter heat bath, which was 

originally defined with respect to the analytically continued Killing vector J^. 
Since t is the proper time of the geodesic observer resting at the origin of 
the static coordinate system (2.6), we conclude that the formula gives the 
probability of a black hole nucleation per Hubble volume and unit time of a 
freely falling observer. 

In order to use the formula (2.15), we should be able to compute the one-loop 
partition functions Z[S 2 x S* 2 ] and ZfS* 4 ]. Now we shall calculate them within 
the path integral approach. 



3 The path integration procedure 



In this section we shall consider the path integral for fluctuations around an 
instanton solution of the Einstein equations R^ u = Ag^ u in the stationary 
phase approximation. We shall largely follow the approach of Gibbons and 
Perry [27]. 



3.1 The second variation of the action 



Our starting point is the action for the gravitational field on a compact Rie- 
mannian manifold Ai, 

M 

whose extrema, 51 = 0, are determined by the equations 

= • (3.2) 

However, if one considers a time evolution along such a family then the action will 
be infinite, which shows that the classical picture does not apply. Even if one uses 
the classical formula for T in this case, one arrives at the quantum result, since 
|a;_|/T=const.~ 1. 
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Let g^ v be an arbitrary solution, and consider small fluctuations around it, 
9fiu —> 9^u + h^ v . The action expands as 

I[g^ + VI = !M + 5 2 I + ... , (3.3) 

where S 2 I is quadratic in h^ u and dots denote the higher order terms. One can 
express 5 2 I directly in terms if h^ v . However, it is convenient to use first the 
standard decomposition of h^ u , 

V = + 4 h 9i*> + V M £, + - - g^aC • (3.4) 

Here <p^ v is the transverse tracefree part, V M 0^ = 0^ = 0, h is the trace, 
and the piece due to £ M is the longitudinal tracefree part. Under the gauge 
transformations (general diffeomorphisms) generated by £ M one has h^ v — > 
h^v + V + V^^ M . The TT-tensor (j)^ is gauge- invariant, while the trace h 
changes as h — > /i + 2V -£ CT . It follows that 

h = h - 2V CT r (3.5) 

is gauge- invariant. For further references we note that £ M can in turn be de- 
composed into its coexact part r/^, for which V M r? M = 0, the exact part V M X) 
and the harmonic piece 

£m = ^ + V m X + £J. (3.6) 

The number of square-integrable harmonic vectors is a topological invariant, 
which is equal to the first Betti number of the manifold Ai. Since the latter 
is zero if M. is simply-connected, which is the case for A > 0, we may safely 
ignore the harmonic contribution in what follows. 

With the decomposition (3.4) the second variation of the action in (3.3) is 
expressed in terms of the gauge-invariant quantities (p^ v and h alone, 

5 2 I = l - (r\ A 2 </v> - 1 (h, A h) . (3.7) 

Here and below we consider the following second order differential operators: 
the operator for the TT-tensor fluctuations 

A 2 </v = -V CT V CT </v - 2iW0 a/3 , (3.8) 
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the vector operator acting on coexact vectors r]^ 

A 1 = -V CT V CT -A, 



(3.9) 



and the scalar operators for h, h, and x 

A = -V a V a , 
A = 3A - 4A , 

A2 = 7 Ao-Ao, (3.10) 

with 7 being a real parameter. Since for A > the manifold M. is compact, 
these operators are (formally) self-adjoint with respect to the scalar product 

(<W^> = / <t>^ v V9d 4 x- (3.11) 

M 



similarly for vectors (77^, 77^) and scalars (x,X)- 

The action S 2 I in (3.7) contains only the gauge-invariant amplitudes (p^ u and h, 
while the dependence on the gauge degrees of freedom £ M cancels. Pure gauge 
modes are thus zero modes of the action. Fixing of the gauge is therefore 
necessary in order to carry out the path integration. To fix the gauge we pass 
from the action S 2 I to the gauge-fixed action 

5 2 I gf = 5 2 I + 5 2 I g , (3.12) 



where, following [27], we choose the gauge- fixing terms as 

5% = 7 / V ff /£ - - p V A VX - i V»h\ , (3.13) 



with 7 and (3 being real parameters. We shall shortly see that it is convenient 
to choose [27] 

,3=-^. (3.14) 



This choice, however, implies that 5 2 I g does not vanish for 7 — > 0. It is often 
convenient to set 7 = 1, in which case (3 — 2. However, we shall not fix the 
value of 7 , since this will provide us with a check of the gauge-invariance of 
our results. 
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Using the decompositions (3.4), (3.6) the gauge-fixing term reads 

6 2 I g = 7 (^, Atf) + JL((h + 2A2x), Mh + 2AM) ■ (3.15) 



Adding this up with 5 2 I in (3.7) one obtains the gauge- fixed action 5 2 I g j. It 
is now convenient to pass from the gauge-invariant variable h defined in (3.5) 
back to the trace h, since with the choice in (3.14) the resulting action then 
becomes diagonal: 



S 2 I 9 f = \ {<T, A 2 </v> + 7 <77 M , A^) (3.16) 
+ \ (x, A A A^ x ) - ±- (h, A h) . 

This action generically has no zero modes, but it depends on the arbitrary 
parameter 7, which reflects the freedom of choice of gauge-fixing. In order to 
cancel this dependency, the compensating ghost term is needed. 



3.2 The mode decomposition of the action 

We wish to calculate the path integral 

Z[g^} = e' 1 J £>[V]£>F P exp (-5 2 I gf ) , (3.17) 

where I = I[g^u\ is the classical action, and the Faddeev-Popov factor is 
obtained from 

1=V FP J D[^\ exp (-8 2 I g ) . (3.18) 



In order to perform the path integration, we introduce the eigenmodes asso- 
ciated with the operators A 2 , Ai and A : 



A 2 tf<g = e fc tf<g, 

A a (p) = A p « (p) . (3.19) 

Throughout this paper we shall denote the eigenvalues and eigenfunctions of 
the tensor operator A 2 by Sk and 0J^, and those for the vector operator Ai by 
a s and rfi*\ respectively. [Later we shall use the symbol s also for the argument 
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of the (^-functions, and this will not lead to any confusion] . Eigenvalues of the 
scalar operator will be denoted by A p , and it will be convenient to split the 
set {\ p } into three subsets, {\ p } = {A ,Aj,A n }, where A = 0, Aj = |A, 
and A n > |A; see Eqs.(3.25)-(3.27) below. Accordingly, the set of the scalar 
eigenfunctions will be split as {a^} = {a^°\a^\a^}. 

Since the manifold is compact, we choose the modes to be orthonormal. This 
allows us to expand all fields in the problem as 

^ = £^, ^ = E^W, (3.20) 

k s 

and 

X = E C ?« (P) » h = Y,C h p a^\ h = J2C~y p) . (3.21) 

V V V 

As a result, the action decomposes into the sum over modes, and the path 
integral reduces to integrals over the Fourier coefficients. 

a) Vector and tensor modes.— Let us consider the mode decomposition 
for the gauge- fixed action in (3.16). This action is the sum of four terms. For 
the first two terms we obtain 

^,A 2 ^) = ±£ efc (<#) 2 , (3.22) 
7^A^)= 7 EW 2 K) 2 ' (3-23) 

s 

These quadratic forms should be positive definite, since otherwise the integrals 
over the coefficients C would be ill-defined. We can see that the quadratic 
form in (3.23) for the vector modes is indeed non-negative definite. Next, the 
expression in (3.22) for the gauge-invariant tensor modes is positive-definite 
if all eigenvalues are positive. If there is a negative eigenvalue, E- < 0, 
as in the case of the S 2 x S 2 instanton background, then it is physically 
significant. The integration over is performed with the complex contour 
rotation, which renders the partition function imaginary thus indicating the 
quasiclassical instability of the system. 

Let us consider now the contribution of the longitudinal vector piece to the 
action (3.16). We obtain 

- ( X , A A Afo> = - £ AAA^ (C*) 2 , (3.24) 

'i *± p 



16 



where \ p = 3\ p — 4A and A^ = 7A P — \ p are the eigenvalues of A and Aq, 
respectively. We note that while X p > 0, the A p and A^ can be negative and 
should therefore be treated carefully. Let us split the scalar modes into three 
groups according to the sign of \ p : 



A a (0) = 0, A = -^A, (3.25) 

A a (t) = ^Aa (i) , \i = 0, (3.26) 

A a {n) = X n a (n \ A n >0. (3.27) 

First we consider the constant mode in (3.25). This exists for any back- 
ground, and for compact manifolds without boundary this is the only nor- 
malizable zero mode of Ao- Since this mode is annihilated by Ao, it does not 
contribute to the sum in (3.24). 

Consider now the scalar modes with the eigenvalue 4A/3 in (3.26). In view 
of the Lichnerowicz-Obata theorem [49], the lowest non-trivial eigenvalue of 
Ao for A > is bounded from below by 4A/3, and the equality is attained if 
only the background is S 4 . Hence the modes in (3.26) exist only for the S A 
instanton, and there can be no modes 'in between' (3.25) and (3.26). In the S A 
case there are five scalar modes with the eigenvalue 4A/3, and their gradients 
are the five conformal Killing vectors that do not correspond to infinitesimal 
isometries. If R^ v = Ag^, a theorem of Yano an Nagano [49] states that such 
vectors exist only in the S 4 case. Let us call these five scalar modes 'conformal 
Killing modes'. Notice that these also do not contribute to the sum in (3.24). 

To recapitulate, the lowest lying modes in the scalar spectrum are the con- 
stant conformal mode in (3.25), which exists for any background, and also 
5 'conformal Killing modes' in (3.26) which exist only for the S* 4 instanton 
and generate the conformal isometries. As we shall see, these 1+5 lowest ly- 
ing modes are physically distinguished, since they are the only scalar modes 
contributing to the partition function. However, they do not enter the sum in 
(3.24). 

For the remaining infinite number of scalar modes in (3.27) (these are labeled 
by n) the eigenvalues A„ and A„ are positive, and it is not difficult to see that 
all the A^'s are also positive, provided that the gauge parameter 7 is positive 
and large enough. To recapitulate, the contribution of the longitudinal vector 
modes to the action is given by 



( X , AoAoA^x) = 7 E ^nXnXl (C£) 2 , (3.28) 
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which is positive definite. We shall see that all modes contributing to this sum 
are unphysical and cancel from the path integral. 

b) Conformal modes.- We now turn to the last term in the gauge-fixed 
action (3.16). Using (3.25)-(3.27) we obtain 

" 1^ < A ' A " H) = I (C « » 2 + - ^ E K (Ci? . (3.29) 

The expression on the right has a finite number of positive terms, correspond- 
ing to the distinguished lowest lying modes, and infinitely many negative ones. 
As a result, an increase in the coefficients makes it arbitrarily large and 
negative, thus rendering the path integral divergent. This represents the well- 
known problem of conformal modes in Euclidean quantum gravity [25]. A 
complete solution of this problem is lacking at present, but the origin of the 
trouble seems to be understood [46]. In brief, the problem is not related to 
any defects of the theory itself, but arises as a result of the bad choice of the 
path integral. If one starts from the fundamental Hamiltonian path integral 
over the physical degrees of freedom of the gravitational field, then one does 
not encounter this problem. The Hamiltonian path integral, however, is non- 
covariant and difficult to work with. One can 'covariantize' it by adding gauge 
degrees of freedom, and this leads to the Euclidean path integral described 
above, up to the important replacement [25] 

h^ih. (3.30) 

The effect of this is to change the overall sign in (3.29), such that the infinite 
number of negative modes become positive. Unfortunately, such a consistent 
derivation of the path integral has only been carried out for weak gravitational 
fields [46] (and for A = 0), since otherwise it is unclear how to choose the 
physical degrees of freedom. Nevertheless, the rule (3.30) is often used also 
in the general case [25], and it leads to the cancellation of the unphysical 
conformal modes. However, some subtle issues can arise. 

For A > the expression in (3.29) contains, apart from infinitely many nega- 
tive terms, also a finite number of positive ones, which are due to the distin- 
guished lowest lying scalar modes. If we apply the rule (3.30) and change the 
overall sign of the scalar mode action, then the negative modes will become 
positive, but the positive ones will become negative. As a result, the path 
integral will still be divergent. It was therefore suggested by Hawking that the 
contour for these extra negative modes should be rotated back, the partition 
function then acquiring the factor , where M is the number of such modes 
[32]. As we know, Af = 6 for the S* 4 instanton, and Af = 1 for any other 
solution of = Ag^ u with A > 0. 
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Unfortunately, this prescription to rotate the contour twice leads in some cases 
to physically meaningless results; the examples will be given in a moment. We 
suggest therefore a slightly different scheme: not to touch the positive modes 
in (3.29) at all and to change the sign only for the negative modes. The whole 
expression then becomes 

- (h - Koh) = i (c5)2 + w ? (c?)2 (C '' )2 ' m) 

We make no attempt to rigorously justify such a rule. We note, however, that 
it is essentially equivalent to the standard recipe (3.30) - up to a finite number 
of modes which we handle differently as compared to Hawking's prescription. 
We shall now comment on this difference. 

When compared to Hawking's recipe [32], the ultimate effect of our prescrip- 
tion is to remove the factor from the partition function. We are unaware 
of any examples where it would be necessary to insist on this factor being 
present in the final result. On the contrary, the examples are in favour of the 
factor being absent. For the S" 4 instanton one has M = 6, such that = — 1, 
and this would render the partition function for hot gravitons in a de Sit- 
ter universe negative, which would be physically meaningless. Next, for the 
S 2 x S* 2 instanton, which already has one negative mode in the spin-2 sector, 
one has Af = 1. As a result, the factor would make the partition function 
real instead of being imaginary, and there would be no black hole pair creation 
I 

These arguments suggest that Hawking's rule should be somehow modified, 
and we therefore put forward the prescription (3.31). Let us also note that our 
rule leads to gauge invariant results - the dependence on the gauge parameter 
7 cancels after the integration. Finally we note that the lowest lying scalar 
modes are physically distinguished, and since they are positive, they should 
be treated similarly to the physical tensor modes. 

To recapitulate, the mode expansion of the gauge-fixed action 5 2 I g f is given 
by the sum of (3.22), (3.23), (3.28), and (3.31): 

5%f = \ E** (Cf) 2 + 7 E(*.) 2 (<7?) 2 + 7 E ^n\nK (C*) 2 

Z k s ^ n 

+ \ (^) 2 + ^ E(^f) 2 + ^ E % (C h n? ■ (3-32) 

In a similar way we obtain the following mode expansion for the gauge-fixing 
term 5 2 I g in (3.15): 
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S% = 7 E(- S ) 2 (^) 2 + A 3 £ (c? - ^C?) 

+ ^^^(2^^ + ^ . (3.33) 
This expression is non-negative definite. 

3.3 The path integration measure 

In order to compute the path integrals in (3. 17), (3. 18) we still need to de- 
fine the path-integration measure. The perturbative measure is defined as the 
square root of the determinant of the metric on the function space of fluctua- 
tions: 

£>[V] ~ yjBet((dh„ v ,dh^)) , D[^\ ~ ^Det((d^, d^}) . (3.34) 

Here it is assumed that the fluctuations are Fourier-expanded and the differen- 
tials refer to the Fourier coefficients, while the meaning of the proportionality 
sign will become clear shortly. Let us first consider -D[^]. It follows from 
(3.4), (3.6) that 

(V, hn = (</v, cjT) + Ai^) + (x, AqAox) + \ (h, h) , 

(Ce) = (^,^) + (X,A X). (3.35) 

Expanding the fields on the right according to (3. 20), (3. 21) and differentiating 
with respect to the Fourier coefficients we obtain the metric for the vector 
fluctuations 



de) = <bf) + (d x , A d x ) = E(^') 2 + E \ idC*Y , (3.36) 

s p 

which yields 

V/Det^,^)) = (j[dcij (jl J\ p dC*J . (3.37) 

Here the prime indicates that terms with \ p = do not contribute to the 
sum in (3.36), and should therefore be omitted in the product in (3.37). To 
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obtain the measure -D[£ M ] we endow each term in the products in (3.37) with 
the weight factor /i 2 / ' \pK\ 

m - (n$ <«?) (n^ /f (n^ yW) . 

Such a normalization implies that 

1= / exp <^)) . (3.39) 

Here fi Q is a parameter with the dimension of an inverse length. In a similar 
way we obtain the measure D[hfj, v ], which is normalized as 




(3.40) 



we shall shortly comment on the relative normalization of _D[/i Mi ,] and -D[£ M ]. 
The result is 




Here the prime indicates that the zero modes of the vector fluctuation oper- 
ator do not contribute to the product. Notice, however, that these modes do 
contribute to the measure D[^]. 

The following remarks are in order. We use units where all fields and pa- 
rameters have dimensions of different powers of a length scale I. One has 
[1/G] = [A] = [// 2 ] = [l~ 2 ]. Eigenvalues of all fluctuation operators have the di- 
mension [l~ 2 ]. The coordinates are dimensionless, while [g^] = [h^] = [I 2 ]. 
For the vectors, [r]^] = = [Z 2 ], and for the scalars [h] = and [x] = [I 2 ]. 
We assume that the scalar, vector and tensor eigenfunctions in (3.19) are 
orthonormal with respect to the scalar product in (3.11). As a result, the di- 
mensions of the eigenfunctions are [<f>$\ = [I], [r]^] — [ct^] = [I -1 ], which 
gives for the Fourier coefficients in (3.20), (3.21) [C*] = [C h ] = [I], [C r <] = [I 2 ], 
and [C*] = [I 3 ]. 

The normalization of D[h^ v ] can be arbitrary, which is reflected in the pres- 
ence of the arbitrary parameter ji in the above formulas. However, the relative 
normalization of D[h^ u ] and -D[£ M ], which is defined by Eqs.(3.39) and (3.40) 
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is fixed by gauge invariance. Had we chosen instead a different relative nor- 
malization, say dividing each mode in (3.38) by 2, then the path integral 
would acquire a factor of 2^0, where Nq is the 'number of eigenvalues' of 
the non-gauge-invariant operator Aq. [The issue of relative normalization of 
the fluctuation and Faddeev-Popov determinants seldom arises, since in most 
cases the absolute value of the path integral is irrelevant]. 

3.4 Computation of the path integral 

Now we are ready to compute the path integrals in (3. 17), (3. 18). Let us il- 
lustrate the procedure on the example of Eq.(3.18), which determines the 
Faddeev-Popov factor V FP . Using S 2 I g from Eq.(3.33) and the measure -D[£ M ] 
from (3.38) we obtain 

(Vfp)- 1 = 11/ dCS exp (- 7 (<7.) 2 (C?) 2 ) (3.42) 

,n/|^*?«p(-^i: k ( 2K q + d) 2 ) , 

which gives 

3.4.I Zero modes 

The factor Qi in (3.43) arises due to the gauge zero modes, for which a s = 
(T = and the integral is non-Gaussian: 




(3.44) 



with the product taken over all such modes. The existence of zero modes of 
the Faddeev-Popov operator indicates that the gauge is not completely fixed. 
This can be related to the global aspects of gauge fixing procedure known 
as the Gribov ambiguity. However, Gribov's problem is usually not the issue 
in the perturbative calculations, where zero modes arise rather due to back- 
ground symmetries. This will be the case in our analysis below. Specifically, 
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the isometries of the background manifold A4 form a subgroup Ti of the full 
diffeomorphism group. Sometimes TL is called the stability group; for the S A 
and S 2 x S 2 backgrounds TL is SO (5) and SO (3) x SO (3), respectively. Since the 
isometries do not change h^ v (in the linearized approximation), their genera- 
tors, which are the Killing vectors K^, are zero modes of the Faddeev- Popov 
operator. 

We therefore conclude that the integration in (3.44) is actually performed 
over the stability group TL. Since the latter is compact in the cases under 
consideration, the integral is finite. In order to actually compute the integral, 
some further analysis is necessary, in which we shall adopt the approach of 
Osborn [44]. First of all, let us recall that all eigenmodes in our analysis have 
unit norm. If we now rescale the zero modes such that the Killing vectors 
Kj = Kj-^ become dimensionless (remember that the coordinates x^ are 
also dimensionless), then the expression in Eq.(3.44) reads 



ni = /ll 



u 2 

^° "Kj\\dCj, (3.45) 



7T 



where now [||ifj||] = [I 2 ] and [Cj] = [I ]. For small values of the parameters Cj 
they can be regarded as coordinates on the group manifold TL in the vicinity 
of the unit element. Since TL acts on M via x^ — > x^(Cj), one has Kj = 
o|r = f^-gfir- However, strictly speaking Cj are not coordinates on the group 
manifold TL but rather on its tangent space at the group unity, such that 
their range is infinite. We wish to restrict the range of Cj, and for this we 
should integrate not over the tangent space but over TL itself. In other words, 
to render the integral in (3.45) convergent we must treat the zero modes non- 
perturbatively, and for this we should replace the perturbative measure Yij dCj 
by a non-perturbative one, djj,{C). 

In general it is a difficult issue to construct the non-perturbative path in- 
tegration measure. However, in the zero mode sector this can be done. We 
note that the measure should be invariant under the group multiplications, 
dfi(CC) = dfi(C), and this uniquely requires that d/j,(C) should be the Haar 
measure for TL. The normalization is fixed by the requirement that for Cj — > 
the perturbative result (3.45) is reproduced. This unambiguously gives 

dfi(C) , (3.46) 

where is computed at Cj = and dfi(C) is the Haar measure of the 
isometry group Ti normalized such that d^iC) — > Yij dCj as Cj — > 0. 



fix 
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3.4-2 The path integral 



Now, using (3.32) and the measure (3.41), we compute the path integral in 
(3.17) - first without the Faddeev- Popov factor V FP : 




(3.47) 



Here the primes indicate that zero and negative modes should be omitted from 
the products. Zero vector modes do not contribute since they are not present in 
the path-integration measure (3.41), and we assume that there are no negative 
vector modes, since otherwise the metric on the space of fluctuations would 
not be positive definite. For tensor fluctuations negative and zero modes are 
present in the measure (3.41), and their overall contribution is collected in 
the factor Q 2 i n (3.47). Let us further assume that there are no zero tensor 
modes, which is the case for the manifolds of interest. If negative modes are 
also absent then Q 2 — 1- If there is one negative tensor mode with eigenvalue 
E- < 0, then 

fi 2 = J=/^exp(-iM^) 2 ). (3-48) 



The integral is computed via the deformation of the contour to the complex 
plane, which gives the purely imaginary result 

Q 2 = , (3.49) 

2i ' 



with the factor of 1/2 arising in the course of the analytic continuation [12]. 

Both the Faddeev-Popov factor in (3.43) and the path integral in (3.47) depend 
on the gauge parameter 7. However, the 7-dependence exactly cancels in their 
product, which provides a very good consistency check. In particular, the 
relative normalization of the integration measures fixed by Eqs.(3.39) and 
(3.40) is important. If we had divided each factor in the mode products in 
(3.38) by a ^ 1, then the resulting path integral would be proportional to 
(Iln a ) ~ a ^ with ( being the C-function of the 7-dependent operator Aq. 
Thus, unless a = 1, the result would be gauge-dependent. 
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We therefore finally obtain the following expression for the path integral in 
(3.17): 

Here Q 2 is the contribution of the negative tensor mode, and Qi is the isometry 
factor. As we expected, the contribution of all unphysical scalar modes has 
canceled from the result. The only scalar modes which do contribute are the 
several lowest lying modes which seem to be physically distinguished. These 
are the constant conformal mode giving rise to the factor //„/ \/2A, and the 
5 'conformal Killing scalars' which exist only in the S A case and give rise to 
the product over %. The next two factors in (3.50) is the contribution of the 
transversal vector modes and the TT-tensor modes. Finally, / = I[g^ v ] is the 
classical action. 

In order to apply the above formula for Z[g^ u ] we need the eigenvalues of 
the fluctuation operators. Now we shall determine the latter for the manifolds 
S 2 x S 2 and S 4 . 



4 Spectra of fluctuation operators 

In this section we derive the spectra of small fluctuations around the S 2 x S 2 
and S* 4 instantons. In the S 2 x S 2 case the problem is tackled via solving 
the differential equations. It turns out that in a suitable basis the system of 
10 coupled equations for the gravity fluctuations splits into 10 independent 
equations. The latter are solved in terms of spin-weighted spherical harmonics. 
In the S A case the equations do not decouple and the direct approach is not so 
transparent. However, the problem can be conveniently analyzed with group 
theoretic methods, which was done some time ago by Gibbons and Perry [27]. 
We shall describe the group theory approach in some detail - with the same 
principal result as in [27]. 



4-1 Fluctuations around the S 2 x S 2 instanton 

Let us consider the metric of the S 2 x S 2 instanton background, 

ds 2 = I ( (dA) 2 + sin 2 ^ (d^) 2 + (d$ 2 ) 2 + sin 2 tf 2 (c^ 2 ) 2 ) • (4.1) 
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It is convenient to set A = 1 for the time being; at the end of calculations the 
A-dependence is restored by multiplying all eigenvalues with A. We introduce 
the complex tetrad 



e 3 = -if d-&2 + -r^rd(p 2 ) , e 4 =4^(^2- 



sini?i 

i 



d<pij , 
■ dip 2 ) ■ 



y/2\ sin ^ 2 / ' V2 \ sin ^ 2 / 

The metric in (4.1) splits as g^ u = e a fi e b v r] a b, where the tetrad metric is 



(4.2) 



10 
10 
1 
10 



(4.3) 



4-1.1 Tensor modes 

First we consider the eigenvalue problem 



(4.4) 



where 



o, 



(4.5) 



We expand <p^ u with respect to the complex basis (4.2), 



(4.6) 



insert this into (4.4) and project the resulting equation onto the basis (4.2) 
again. Remarkably, the system of 10 coupled equations splits then into 10 in- 
dependent equations for the 10 tetrad projections <3> a &. A partial explanation 
of this fact is the existence of two different parity symmetries acting indepen- 
dently on the two spheres. 



It is convenient to introduce the operator 

-Ar a i 9 2 n d n . cottf 1 

V[s, (p] = — + cot ■& — + 2is — — r + 



d 2 



d$ 2 



sin$ sin i? d(jy 



-s 2 cottf, (4.7) 
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whose eigenfunctions are the spin- weighted spherical harmonics s Yj m [29], 

V[s, tp] s Y jm (#, 0) = (s 2 - j(j + 1)) s Y jm (&, 0). (4.8) 



Here j and m are such that j > |s| and —j<m<j. One has B Yj m = 
for j < |s|. [Notice that we use the bold-faced s for the spin weight.] The 
following relations between harmonics with different values of the spin weight 
s are useful: 



£±[8, 0, <p] s Y jm = ±J(j ± s) (j =f s + 1) STl Y jm , (4.9) 

where 

M = ^ T |r ± s cot . (4.10) 
cw smw a0 

The harmonics for a fixed s form an ortho normal set on S 2 . 
Using the above definitions, Eqs. (4.4) can be represented as 

(f>[s^,^ 1 ,v 5 i]+P[s 2 b ,^,^ 2 ]-(s^) 2 -(s 2 6 ) 2 + 2 + £ ) $ ab = 0, (4.11) 

where 1 < a, b < 4 (no summation over a, 6). Here the nonzero elements of 
the symmetric matrices s^ b and s 2 fe are 

s n = — s 22 — 2, s 13 = s 14 = — s 23 = — s 24 = 1 , 

s 2 3 = — s 2 4 = 2, s 2 3 = s 23 = — s 2 4 = — s 24 = 1 . (4-12) 
The solution to Eqs. (4.11) reads 

§ab = C ab sl b Yj imi (A, fi) B^Y j2ma (^ 2 , ¥2), (4.13) 

with C a b being integration constants. The eigenvalue, e, is the same for all 

e = j 1 (j 1 + l)+h(j 2 + l)-2. (4.14) 

This is essentially the sum of squares of the two SO (3) angular momentum 
operators acting independently on the two spheres. 

Let us now count the degeneracy of the modes. For this one should take 
into account the additional conditions (4.5), which gives algebraic constraints 
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for the coefficients C a b- We consider first the trace condition 0^ = 0. In the 
language of the tetrad projections this reduces to $ 12 + $34 = 0, or equivalently 
to 

Ci2 + C 34 = 0. (4.15) 

Hence only 9 out of the 10 constants C ab are independent. 
Next, we consider the Lorenz condition V CT 0^ = 0. This implies 

t [S^, 0i, tpi] $„l + C + [sl 2 , 0!, Vl ] $ a2 
+ t [S^ 3 , 2 , </? 2 ] $ a3 + C + [s 2 a4 , 2 , <^ 2 ] $a4 = (4.16) 

(no summation over a). Inserting the solution (4.13) and using the recurrence 
relations in (4.9), these conditions reduce to algebraic constraints 

Ki Cn - OL\Cv2 + Ct 2 (Ci3 - C U ) = , 

ai C12 - K1C22 + «2 (C23 - C24) = , 

«1 (C13 - C23) + «2 C33 - «2 C34 = , 

«i (C14 - C 24 ) + « 2 <^34 - k 2 C 44 = . (4. 17) 

Here a L = ^j L (j t + l) (with 1 = 1,2) and « t = yj(j t + 2)(J t - 1) for j t > 1 
while K t = for j t = 0. 

For j L > 2 (which corresponds to quadrupole or higher deformations of each 
sphere) none of the coefficients a L , k l vanish, and the algebraic constraints 
(4.17) reduce the number of independent coefficients C ab to five. This gives 
the degeneracy d: 

ji > 2, 32 > 2, d = 5 (2j4 + l)(2j 2 + 1) . (4.18) 

The situation is different for small values of j L . Consider, for example, the 
ji — 32 — sector. Since the harmonics s Yj m vanish for 3 < |s|, we must set 
in the solution (4.13) all C a b's to zero, apart from C\ 2 = —6*34. The Lorenz 
condition (4.16) is then fulfilled. As a result, there is only one independent 
integration constant, which yields 

ji=32 = 0, d=l. (4.19) 
Notice that in this case e = —2. 

In a similar way one can consider the sector where j\ = and j 2 = 1 (or 
ji = 1 and 32 = 0), in which case e — 0. One discovers then that the Lorenz 
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constraints (4.16) require that all non-trivial coefficients C a b mush vanish. As 
a result, 

j! = 0, j 2 = 1, or j 2 = 1, j! = 0, d = 0, (4.20) 



which shows that there are no zero modes. 
Next, 

h =32 = 1, d = (2 3l + l)(2j 2 + 1) = 9; (4.21) 

and finally 

ji > 2, j 2 = 0, d = 2 3 \ + 1; 

3 \ > 2, j 2 = 1, d = 9(2j! + 1), (4.22) 

which conditions are symmetric under interchanging ji and j 2 - 

To recapitulate, the spectrum of the tensor fluctuations contains one negative 
mode and no zero modes. 

4-1-2 Vector modes 

Let us now consider the eigenvalue problem 

(-V a V Q -A)^ = a^ (4.23) 

subject to the condition 

V M ?f = (4.24) 

for the vector fluctuations on the S 2 x S 2 background. We again expand the 
fluctuations with respect to the basis (4.2), 



V, = e a ^ a , (4.25) 



insert this into (4.23), and project back to the tetrad. Similarly as in the tensor 
case, the equations decouple to give 

(V[sl A, V?i] + V[s 2 a , 2 , V2] + 1 + *) = , (4.26) 
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where 1 < a < 4 (no summation over a), and nonzero coefficients read s\ = 
—si = s§ = — s| = 1. The solution is 

= C a si*Wtfl^l)s^ 2 m 2 (#2,¥> 2 ), (4.27) 

with C a being integration constants, and the eigenvalue is the same for all \l/ a : 

<j = j 1 (j 1 + l)+j 2 (j 2 + l)-2. (4.28) 

The Lorenz condition, V ' a rf = 0, reads 

£- [si, v?i] + £ + [s^, v?i] ^2 
+ 1 [si, 2 , V? 2 ] ^3 + ^ + [s^, ^2, <^ 2 ] ^4 = 0, (4.29) 

which reduces upon inserting (4.27) to the algebraic condition 

VJi(Ji + 1) (Ci - C 2 ) + VSJ2TT) (C 3 - C 4 ) = 0. (4.30) 

This allows one to count the degeneracies: 

ji > 1, 32 > 1 (o- > 0), d = 3 (2ji + l)(2j 2 + 1); (4.31) 

and also 

Ji>2, j 2 = (<r>0), d^j'iO'i + l); 
Ji = l, J 2 = (a = 0), d = 3; 

Ji = 0, j 2 = (a = -2), d = 0. (4.32) 

These results are symmetric under ji <-> j'2, hence there are no negative modes, 
there are six zero modes corresponding to the six Killing vectors of S 2 x S 2 , 
and the rest of the spectrum is positive. 

4-1-3 Scalar modes and the orthogonality conditions 
The eigenvalue problem for the scalar modes, 

-V a V a h = Xh, (4.33) 

reduces to the equation 

(V[0, A, <pi] + V[0, 2 , <p 2 ] + A) h = , (4.34) 
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whose solutions are 

h = >W (^1^1)^2(^2,^2) 



(4.35) 



(for s = the spin-weighted spherical harmonics coincide with the usual spher- 
ical harmonics). The eigenvalues are just 

A = Ji(ji + l)+J 2 (j 2 + l), (4.36) 

and the degeneracies are 

ji > 0, 32 > 0, d = (2j! + l)(2j 2 + 1). (4.37) 

We have obtained the spectra of all relevant fluctuation operators. Although 
the eigenfunctions are complex, one can pick up their real part in a way that 
is consistent with the orthogonality conditions. For example, for the scalar 
modes one considers 

®(h) = ^ Y jimi Y nm2 + c.c , (4.38) 



and one can easily see that the modes 9ft (/i) with different quantum numbers 
{j\m\3\m-i) are orthogonal with respect to the scalar product defined in Eq. 
(3.11). 

For the vector modes ty a the procedure is slightly more complicated, since 
the tetrad metric i] ab is not diagonal. In addition, harmonics s Yjm for differ- 
ent values of the spin weight are not orthogonal. Consider, however, the real 
combinations 

= ±±1 e a ^ a + c.c , (4.39) 



where \l/ a has quantum numbers (jimij^r?^) . Consider 3t(r)(p) and 3t(r)^) 
with different quantum numbers. Their scalar product (defined in Eq. (3.11)) 
can be computed using the relations 

Vab = efeS g, u , 5 ab = e/(e^)* g, v , (4.40) 

which gives 

(K«)),K(^)) = E(^ 1) ^i 2) *) (4-41) 

a 

+i(*S 1) ,^ 2) )+i(*? ) ,*i 2) ) + c.c. 
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Table 1 

Spectra of fluctuations around the S 2 x S 2 instanton 



operator 


eigenvalue 


degeneracy 




A 2 


-2A 
2A 

+ 1) - 2) A 

j(j + 1)A 
(Ji(ii + 1)+J2(i2 + 1)-2)A 


1 

9 

2(2j + 1) 
18(2.7 + 1) 
5(2ji + l)(2j 2 + 1) 


J > 2 
J>2 
Jl, J2 > 2 


Ai 




+ 1) - 2)A 
(Ji(ii + 1)+J2(i2 + 1)-2)A 


6 

2(2j + 1) 
3(2ji + l)(2j 2 + 1) 


J>2 

J1,J2 > 1 


A 


(ji(j'i + l)+i 2 (j2 + l))A 


(2ji + l)(2j 2 + 1) 


3i, h > 



Here each term in the sum J2 a (^ii\ ^^*) * s a bilinear combination of spin- 
weighted harmonics with the same value of the spin weight, such that the or- 
thogonality relation holds. Next, integrating by parts and using the recurrence 
relations (4.9) one can show that the remaining term in the scalar product, 
i(^f^\ +^(^3 1 \ ^i 2 '*) + c -c 5 vanishes. This shows that vector modes ^(r/^) 
with different quantum numbers are orthogonal. 

A similar procedure can be carried out for the tensor modes. Hence for all 
eigenmodes considered above one can choose the real part in such a way that 
the orthogonality condition holds. This is a manifestation of the fact that the 
fluctuation operators are self-adjoint. We finally restore the dependence on A 
and summarize the results of this section in Tab.l. There is one negative mode 
in the spectrum, and this plays a crucial role in our analysis. The corresponding 
deformation increases the radius of one of the spheres, shrinking at the same 
time the second one. 



4-2 Fluctuations around the S 4 instanton 

The S 14 instanton can be viewed as the four-dimensional sphere with radius 
y // 3/A in five-dimensional Euclidean space E 5 . Although the corresponding 
eigenvalue problem for fluctuations was considered in [27], we have reanalyzed 
it for the sake of completeness (with the same result) and shall present below 
the key steps of our analysis. The problem essentially reduces to studying 
representations of SO(5) [31,1,2,8], whose Casimir operator can be related 
to the invariant Laplacians on S* 4 with the help of the projection formalism 
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[34]. We shall therefore first outline the group theory part by summarizing 
the relevant facts about representations of 5*0(5). We shall work on the unit 
4-sphere rescaling at the end the eigenvalues by the factor A/3. 



4-2.1 Representations of SO (5) 

The unit sphere S 4 in E 5 is defined in Cartesian coordinates by the equation 
Yf a=1 (x a ) 2 = 1. It is convenient to use the complex coordinates £ ±x = (x 1 ± 
ix 2 )/\^2, £ ±2 = (x 3 ±ix 4 ) I y/2, £° = x 5 . We shall not distinguish between lower 
and upper indices, = £\ In these new coordinates the defining quadratic 
form reads Z)i=-2'C l ^~* = 1> which is annihilated by the generators of S0(5): 

whose commutation relations are 

[Yj, Yf] = StY? - S\Y} + 5~ l Yi k - S'UYT* = C?j kl Y? . (4.43) 

Since Y- = —Y/, the independent generators can be chosen to be those for 
—% < j. Y] 1 and F 2 2 generate the Cartan subalgebra, while YJ and Y? for 
—i<j<i are the raising and lowering operators, respectively. One has 

[Y7, Yf] = at(i)Yf, (4.44) 



where 

of (0 = " *i + £-1 " *l* ( 4 - 4 5) 



determine the root vectors with the components af = (ctf (1), ctf (2)). The 
roots corresponding to the four raising operators are a\ = (—1,1), a 2 , = (0, 1), 
a 2 _ l = (1, 1), and = (1, 0). 

Irreducible representations of SO (5) are characterized by two numbers denoted 
by m = (mi,m2), where r«2 > mi and both m\ and m 2 are either integer 
or half-integer. The highest weight vector %p m is annihilated by all raising 
operators, YJip m = for i > j > —i, and it is an eigenvector of the Cartan 
subalgebra generators, Yfy m = miifj m , i — 1,2. Using these properties and 
also [Y^Y/] = Yl — Yj — 25 t _jYl j , one finds the eigenvalues of the Casimir 
operator, 

C^ m = W Y J Y ^rn = C m ifj m , (4.46) 
i,3 
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where 



C m = mi (mi + 1) + m 2 {m 2 + 3) (4.47) 

is the same for all vectors of the representation. The dimension of representa- 
tions is given by 

dim(ra) = JJ(o;, r + m) / r) . (4.48) 

a a 

Here the product is over the four root vectors described above, and r = 
\ Ha a — {\i f )■ One nas r + m = +mi, § +m 2 ), and ( , ) is the scalar with 
respect to the Cartan metric = —C^ kl C^ pq = 65ij (here i,j = 1,2). As a 
result, 

dim(m) = ^ (2m x + l)(2m 2 + 3)(m 2 - m 1 + 1)(toi + m 2 + 2) . (4.49) 
6 

4-2.2 Scalar modes 

Using Eqs.(4.47),(4.49) one can find the spectra of the relevant fluctuation 
operators. It is now convenient to pass back to the Cartesian coordinates 
x a = x a {a = 1, . . . 5), such that the sphere S 4 is determined by the condition 
r = ^Jx a x a = 1. Let n a = x a /r be the unit normal to the sphere. The (anti- 
hermitean) generators of SO (5) in Cartesian coordinates are given by M a t, = 
n a db — n b d a , and the Casimir operator is C = —\(M ab ) 2 = — \ Y.ab^Mab) 2 ■ 

Let us define the projection operator = 5 a b — njib = P ab = P fe a , which can 
be thought of as the induced metric on the sphere. In what follows we shall use 
the projection formalism [34] to describe geometrical 4-objects tangent to the 
sphere in terms of 5-objects of the embedding space. For example, a 4-vector 
r]^ can be described as a 5- vector r\ a subject to the condition n a i] a = 0. The 
covariant derivative of a tensor is obtained by taking the partial derivative 
and then projecting all the indices down to the sphere. For example, V a i]b — 
(d p r] q )P^P b q . One has n a = n a , while for objects tangent to the sphere 5-indices 
can be raised and lowered either with P ab or with 5 a b- The curvature tensor is 
given by R psqt = P pq P st - P pt P sq - 

Consider first scalar fluctuations. The covariant Laplacian for a scalar field h 
can be expressed in terms of the angular momentum operator as 

Uh = P ab d a (P b c d c h) = ^(M ab ) 2 h = -Ch. (4.50) 
Scalars transform according to the (0, j) representations, which correspond to 
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the Young tableaux |i |a | . . . | il and can be represented in terms of homogeneous 
polynomials on S* 4 as 

/' /'.-, • (4.51) 



Hence, the eigenvalues of the Casimir operator in Eq. (4.50) are C( 0j ) = 
j(j + 3), which gives the spectrum of the scalar eigenvalue problem, A h = \h 
with A = -f □: 

A = ^j(j + 3), d=^(2j + 3)(j + 2)0' + l), J>0. (4.52) 



4-2.3 Vector modes 

Consider a tangent vector field i] s = Pgf] a . The invariant Laplacian reads 

n Vs = P ab d a (P b c d cVp PV)P« = ^(M ab ) 2 Vs + 2(d aV a )n s + Vs . (4.53) 

Here the last two terms on the right can be related to the contribution of 
the spin operator. Under general SO(5) rotations a vector rj(x) transforms 
into fj(x) = R?](R~ 1 x), where R = exp(u ab S ab ) with uj ab = —oo ba being the 
rotation parameters and S ab = (S a b) pg = S^S^— 5 b 5l. For |o; ab | 1 one obtains 
77 — 77 = u ab (M ab + S a b)i], such that S ab can be identified with the spin operator: 
(-S' a6 77) s = (S ab ) s p r] p . As a result, 

n Vs = Q(M a6 + S ab ) 2 + 3^)% = (-C + 3) Vs , (4.54) 

where the Casimir operator is now the square of the total angular momentum. 
The general vector harmonics on S* 4 can be obtained by considering the prod- 
uct of a constant vector in E 5 with scalar harmonics on S* 4 . Such a product 
decomposes into three irreducible pieces, (0, 1) <g> (0, j) = (1, j) © (0, j + 1) © 
(0, j — 1), which can be visualized as 



\2 1. rrrj] = nnrr. . . i j i © Ri i. . . i j i © fi n i. . . im| . (4.55) 



The first term on the right here is the (1, j)-piece, and in the language of 
homogeneous polynomials it reads 

Vs = V[sa](a 1 ...a,. 1 )n a n ai . . . n a ^n a ^ , (4.56) 



where 77[ sa ]( ai ... a _ 1 ) is traceless with respect to any pair of indices. This is 
manifestly tangential and coexact. As a result, the eigenvalues of the Casimir 
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operator are C^j) = j(j + 3) + 2, and the spectrum of the vector eigenvalue 
problem Ai?y s = (— — A)r/ S = ar/ s in the sector where d a rf = n a r] a = is 
given by 

<r=!0'0' + 3)-4), d=^j(j + 3)(2j + 3), j>l. (4.57) 

One can also directly verify that the harmonic rj s in Eq.(4.56) fulfills the 
condition j(M a b) 2 r] s = —j(j + 3)i] s . It follows then from Eq. (4.53) that Or) s = 
— + 3) — l)i] s , and this again yields the spectrum in Eq. (4.57). The 
correct degeneracy can be obtained by counting the independent components 

Of J7[«a](oi-ai-i) W 

The remaining two pieces in (4.55), when represented in terms of the polyno- 
mials on S* 4 , can be related to the exact tangential and the normal components 
of the vector field. 



4-2.4 Tensor modes 

For a symmetric tensor h pq = PpP b h a b a direct calculation gives 

Uh m + 2R psqt h st = P ab d a {P c b {d^ m )Pl_Pl)PM + 2(P pq P st - P pt P sq )h st 
= X -(M ah fh pq + 2n p (d a h aq ) + 2n q (d a h ap ) + 25 pq h a a 
= {^-{M ab + S a6 ) 2 + 6) h pq ee (-C + 6) h pq . (4.58) 

Here the spin operator is defined in the same way as for vectors, which gives 
(S a b/i) pg = (S a b) ph sq + (S a b) q h sp . The general tensor harmonics on S* 4 are 
obtained by the direct products (0, 2) <g) (0, j) = (0, j + 2) © (1, j + 1) © (0, j) © 
(2, j)©(0, j + l)©(l, j — 1)©(0, j — 2). Again this can be visualized by Young's 
diagrams and represented in the language of homogeneous polynomials. The 
transverse and tracefree harmonics tangent to the sphere correspond to the 
(2,j) piece, whose explicit representation is 

h pg = h [pa][qbKai ...a J - 2 )n a n b n ai . . . n a ^ 2 . (4.59) 

Here h\p a ^ q b]( ai ...a _ 2 ) is traceless with respect to any pair of indices and is sym- 
metric under interchange of the [pa] and [qb] pairs. As a result, the eigenval- 
ues of the Casimir operator are C( 2 j) = j(j + 3) + 6. This gives the spec- 
trum of the tensor eigenvalue problem A 2 h pq = eh pq in the sector where 
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Table 2 

Spectra of fluctuations around the S 4 instanton 



operator 


eigenvalue 


degeneracy 




A 2 


fi(i + 3) 


|(j-l)(j + 4)(2j + 3) 


i>2 


Ai 


#(i0' + 3)-4) 


5J'0' + 3)(2j + 3) 


i > i 


A 


fi(i + 3) 


i(j + l)(i + 2)(2j + 3) 


3 > 



<9 a /i afe = n a /i ab = = 0: 

e = |j'0' + 3), rf=^(j-l)(j+4)(2j + 3), j>2. (4.60) 

The same result can be obtained by directly verifying that h pq in Eq. (4.59) 
fulfills the condition \(M ah ) 2 h pq = + 3)h pq . 

The other tensor harmonics in the expansion of (0,2) ® (0, j) correspond to 
the exact and coexact pieces of the longitudinal vector part of the 4-metric, 
to those of the 4- vector h 5fl , and to the trace [27]. 

We summarize the results of our analysis in this section in Tab. 2. Notice that 
the scalar and tensor eigenvalues are the same (for j > 2), while the vector 
spectrum is shifted by a constant. 



5 Partition function 

Now we are able to derive the explicit expressions for the one-loop partition 
functions for fluctuations around the S 2 x S 2 and S 4 instantons. The corre- 
sponding formula was obtained in Eq.(3.50) above. It is convenient to pass 
from /i to the dimensionless normalization parameter /i via the rescaling 

/i o = v / A/i - (5.1) 

The one-loop partition function for gravity fluctuations around an Euclidean 
background then reads 

7\ l ^° n fi2 / Det/A i -/ (r o\ 

Here the first two factors on the right are the contributions of the scalar modes. 
The factor /x /\/2 is due to the constant conformal mode, which is always 
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present, and fio is the contribution of the 5 scalar modes with eigenvalue 
4A/3 which exist only for the S 4 instanton (see Tab. 2): 



12 1_ 
3 fx 



(5.3) 



For any background other than S A one has f^o = 1- As was discussed above, 
other scalar modes do not contribute to the partition function. 



The factor Vt 2 in (5.2) is the contribution of the negative tensor mode, 



n, = 



2iJ\eJ 



(5.4) 



For the S 2 x S 2 instanton there is one such mode with e 
S A case all tensor modes are positive and fl 2 — 1- Next, 



d 



Vol{H) , 



—2, while in the 



(5.5) 



is the isometry factor. If the background has no isometries then Qi — 1. 

The determinants in Eq. (5.2) are the contributions of the positive vector and 
tensor modes. One has 

^Det'Ax = =exp(~Ci(0)-^(lB^)Ci(0)) , (5.6) 



where the ^-function for the positive, transverse vector modes is 



(57) 



Similarly for the positive, transverse traceless tensor modes: 

/ Det'A 2 = ( I]' x f^] = exp (-±£(0) - ±(\nfi) C 2 (0)) (5.8) 



with 



c,w = E'(^ ; 



(5.9) 
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The last factor in Eq. (5.2) is the classical contribution, with / being the action 
for the background. Let us now apply these formulas. 



5. 1 The S 2 x S 2 instanton 



The classical action is I[S 2 x S 2 ] = — 2n/AG, and according to Tab.l, 



Consider now the isometry factor Qi in (5.5), which is due to the back- 
ground SO(3)xSO(3) symmetry. Each of the two SO(3) groups can be pa- 
rameterized by matrices Uik = exp(eikjCj). The invariant metric on the SO(3) 
space is = |tr(<9j[/<9fc[/ _1 ) — > 8^ for Cj — > 0. The Haar measure is 
d/i(C) = y /Mg^dC 1 dC 2 dC 3 , and the volume Vol(SO(3))= J dpi(C) = 8vr 2 . 
For later use, we reproduce this result in a different way. The measure for a 
(compact, semi-simple) Lie group Q can be represented as the product of the 
measure for the maximal subgroup 7i and that for the coset G /Ti. This implies 
that 

Vol{G) = Vol{H) x Vol{G/H) . (5.11) 



In particular, Vol(SO(3))=Vol(SO(2))xVol(S 2 ), where Vol(SO(2))= 2tt, and 
the volume of the S 2 coset with unit (due to the normalization of the measure) 
radius is Vol{S 2 ) = 4tt. As a result, Vol (SO (3)) =2% x 4tt = 8vr 2 . 

When acting on S 2 , the SO(3) generators generate rotations in the three 

orthogonal planes of the embedding Euclidean 3-space. Let -^r be the genera- 
tor of rotations in the XY-plane, such that the azimuthal angle of the spherical 
coordinate system changes as ip — > y? + C3. Then the norm 1 1 | is the square 
root of 

S 2 xS 2 



Obviously, the norms 



and 1 1^| 



and those of the generators of the 



second SO (3) factor are the same. Hence, 





A 



7T 



dip 



(Vol(SO(3))Y 



647T 4 (/i ) 12 

27(AG) 3 



(5.13) 
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Consider now the positive modes. The ^-function associated with the positive 
vector modes is (see Tab.l) 

Cl(s) = y 2 ^ + *) + yy 3( 2jl + i)(2 j2 + 1) 



This can be represented as 

Ci(s) = 4 s (2 C(2, -9|s) + 3 Z(l, -10\s)), (5.15) 

where the following two functions have been introduced: 



C(Ms) = E r/n- i"> ( 5 - 16 ) 

1 ' ' ^{(2j + l) 2 + ^} s ' 



oo oo 



z(t, ^ )= ££{(%+i)»+(%+i)'+-}- (5 - 17) 

These functions are studied in detail in the Appendix. Similarly, using the 
results of Tab.l one obtains the ^-function for the positive tensor modes 

C 2 (s) = 9 x 2~ s + 4 s (2 C(2, -9|s) + 18C(2, -l|s) + 5Z(2, — 10|s)). (5.18) 



The following relation implied by the definitions in (5.16), (5.17), will be useful: 
Z(l, -10\s) = Z(2, -10|s) + 6C(2, -l\s) + 9 x 8~ s . 



5.1.1 The scaling behaviour 

Before we proceed further, it is very instructive to pause and check whether 
the expressions above agree with the general formulas for the scaling behaviour 
of effective actions. We shall follow the approach of Christensen and Duff [13], 
who relate this scaling behaviour to 



N = 18Q (4?r)2 / (iW^ pCT + Q3QA 2 )Vgd 4 x , 

Ni = m ^ )2 I (-lliW^* 7 + 984A 2 ) v ^rf 4 x , 

N 2 = ^o^a / (l89R, vpa R^ - 756A 2 )^rf 4 x . (5.19) 

Here N is the 'number of eigenvalues' of the scalar operator A — 2A acting 
on a manifold with R^ u = Ag^. JVi is the number of eigenvalues of the vector 



40 



operator Ai acting in the space of all vectors, that is, including both trans- 
verse and longitudinal fluctuations. Finally, N 2 counts both transverse and 
longitudinal eigenstates of the tensor operator A 2 , with the only requirement 
that the fluctuations must be traceless. 

Let us apply these formulas to the S 2 x S 2 background. The volume of the 
manifold is V S 2 xS 2 = (4vr) 2 /A 2 , while R^ pa R^P a = 8A 2 . As a result, 

161 224 21 

N ° = lf N ' = T- (5 - 20) 



Now let us obtain the same result via a direct evaluation of the (^-functions. 
First we consider the scalar case. Using the results of Tab.l, the operator 
A — 2A has one negative mode, six zero modes, while the rest of the spectrum 
is positive and gives rise to the C-function 

Co(s) = 4 s (2 C(2, -9|s) + Z(l, -10\s)) . (5.21) 



Hence the number of all eigenvalues is 7 + Co(0). In order to compute Co(0), 
we use the results of the Appendix, where the following values are obtained: 

C(M0) = ^-j"-* 2 , (5.22) 
Z(M0) = ^-^ + 2^+(i,-^ + i|. (5.23) 

This gives for the ^-functions in (5.15), (5.18), (5.21) 

Co(0) = -±f, d(0) = -f, C 2 (0) = f. (5.24) 



Using these, the number of scalar eigenvalues is N = 7 — ^ = which 
agrees with (5.20). 

Next, the vector operator Ai has 6 zero modes, such that the number of its 
eigenvalues in the transverse sector is 6 + Ci(0)- Now, one should take into 
account also the longitudinal vectors, which are gradients of scalars. It is not 
difficult to see that if V M x is an eigenvector of A x , such that AiV^x — °"V M x, 
then (A — 2A)x = ax We see that the eigenf unctions of A — 2A are in 
one-to-one correspondence with the longitudinal vectors. The number of the 
latter is therefore N — 1, where the one is subtracted because the ground 
state scalar eigenfunction is constant, which vanishes upon differentiation. We 
therefore conclude that AT 1 =6 + Ci(0)+A" -l = 6- f + ^- l = 2 ^, 
which also agrees with (5.20). 
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Finally, the number of traceless eigenvalues of A2 is 1 + £2(0) (here the one 
is the contribution of the negative mode) plus the number of longitudinal 
traceless tensor harmonics 0^ = V M £„ + V„£ M — \ g^ v V ' p ^ p . 

Now, if = a£ u then for (ffc v associated with £ M one has A 2 0j^, = o^ v . 
Hence, the number of longitudinal tensors is determined by the number of 
vectors, which gives = 1 + C2(0) + (N\ — 6). Here six is subtracted because 
the six Killing vectors do not contribute to the tensor spectrum, since for 
Killing vectors one has 0^ = 0. We therefore obtain jV 2 = l + y + ^-6 = y, 
which again is in perfect agreement with (5.20). 

The overall scale dependence of the partition function is given by the factor 
(Ho) N2+N °- 2Nl . For the S 2 x S 2 instanton one has N 2 + N - 2N ± = and 
we shall shortly see that this agrees with our analysis. 



5.1.2 The partition function Z[S 2 x S 2 ] 

It is now a simple task to insert the formulas above into the expression for the 
partition function. We obtain 

= exp (C(0) + In n 2 C(0)) , (5.25) 



where 

C(s) = l - (( 2 (s) - Ci(s)) = -9 x 2"' + A s Z(2, -10\s) . (5.26) 



Using the values Z(2, — 10|0) = and Z'(2, — 10|0) = T = -18.3118 (see 
Eq.(A.51) in the Appendix) we find 



/Det'Al 1567 M2 

5;^ =2^„«e T . (5.27) 



Finally, taking into account the contributions of the negative, zero, and scalar 
modes computed in (5.10), together with the classical term, we obtain 



Z | S . xSl= _,^M e . (5.28) 
1 J 256 7r 4 /!o° V Det' A 2 v ; 



(—) 

\AGJ 



= -i 0.3667 x (AG)Vo 45 ex P 
This is our final result in the S 2 x S* 2 sector. 
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5.2 The S 4 instanton 

The classical action is I[S 4 ] = — 3n/ AG. Using the results in Tab. 2 we find 

Q ° = I = — I , ^2 = 1- (5.29) 




Let us consider the symmetry factor f^. The isometry group is now 7i=SO(5), 
and this can be represented by matrices = exp(Cjfc), where = —Cki, 
i, k = 1, . . . 5. The 10 generators jSr- generate rotations of S 4 in the 10 orthog- 
onal planes of the embedding Euclidean 5-space. Let be the generator of 
rotations in the XY-plane, such that the standard azimuthal angle changes as 
ip — > ip + C\2- The norm || g§^\ | is the square root of 

. 8 8 . 1 f _ , 4 97T 

< v d^ ) = ^G{ 9w ^ 9d x = m^c ' (5 - 30) 



which applies also to the the norms of the remaining 9 generators. 

The volume of SO (5) can be computed by directly constructing the invari- 
ant metric and the Haar measure with the use of the matrix representation 
Uik = exp(Cjfc). The measure should be normalized such that for — > it 
reduces to lli<fe dC^. However, it is much simpler to use the coset reduction 
formula (5.11). One has SO(5)/SO(4) =1 S 4 and SO(4)/SO(3)=,S 3 , such that 
Vol(SO(5))=Vol(S 4 ) x Vol(S 3 ) x Vol(SO(3)). We know that Vol(SO(3))=8n 2 , 
while the volumes of unit S 3 and S 4 are 2n 2 and 87r 2 /3, respectively. As a re- 
sult, 1/o/(SO(5))=128tt 6 /3. Summarizing, 




d 



dip 



10 



9 \ 5 128tt 6 (/i ) 



20 



Vol(SO(5)) = ( — J — (5.31, 



Let us consider the positive modes. The ^-function associated with the positive 
vector modes is (see Tab. 2) 

C^^sj^f^. (5.32) 



This can be written as 



Ci(s) = i3 s Q(l,-4,0| S ), (5.33) 
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where the following function has been introduced 



Similarly, using the results of Tab. 2, one obtains the ^-function for the positive 
tensor modes 

C 2 (s) = |3 S Q(2,0,-4| S ). (5.35) 

Finally, consider the scalar operator A — 2A. According to Tab. 2, its eigenval- 
ues, measured in units of A, are given by + 3) — 6)/3, and the degeneracy 
is (j + l)(j + 2)(2j + 3)/6 with j > 0. Hence, the (^-function for the positive 
scalar modes is 

Co(s) = ^3 s Q(2,-6,-4|s). (5.36) 



5.2.1 The scaling behaviour 

Let us again check the consistency with the general formulas for the scaling 
behaviour of quantum fields (for fluctuations around S* 4 this was done by 
Christensen and Duff [13]). Applying again the formulas in (5.19), where now 
the volume of the manifold is V s * = 24tt 2 /A 2 , while R^ pa R^ pa = 8A 2 /3, one 
has 



On the other hand, using the result of the Appendix, 



Q(k, u, c|0) = - 1 - k 4 - 2k 3 - (c + ^)k 2 (5.38) 

one obtains for the (^-functions in (5.32), (5.33), (5.35) 

C„(0) = -|, C.(0) = -^, «0) = -5i. (5.39) 



Now, since the spectrum of A — 2A contains six non-positive modes, one has 
N = 6 + Co(0) = 6 - § = 4 g, which agrees with (5.37). Next, A x has 10 
zero modes, such that there are 10 + Ci(0) transverse vector eigenstates, plus 
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(Nq — 1) longitudinal ones (the constant scalar mode gives no contribution). As 
a result, jV 1 = 10-^- + ^- l = ^f, which agrees with (5.37). Finally, there 
are N 2 = C2(0)+A 1 — 15 traceless tensor modes, where 15 is subtracted because 
10 Killing vectors and 5 conformal Killing vectors of S A do not contribute to 
the longitudinal tensor modes. One obtains N 2 = — |j + ^ — 15 = — ^, which 
again agrees with (5.37). 



The overall scale dependence of the partition functions is expected to be 

571 
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M Na+N °- 2Nl , where N 2 + N - 2N t = — 571 



5.2.2 The partition function Z[S 4 } 

Let us now obtain the partition function. One finds 



'Det'Ai 
Det'A 2 



exp(C'(0) + ln/ioC(0)) , (5.40) 



where 



C(s) = \ (Us) - Ci(s)) = 3 s (A Q(2, 0, -4| S ) - \ Q(2, -4, 0| S )) . (5.41) 



One has C(0) = ff and C'(0) = T x = 6.1015 (see Eq.(A.36) in the Appendix). 
This yields 



/Det'Ai sea T , 



Finally, collecting the contributions of the negative, zero, and scalar modes 
computed in (5.10), together with the classical term, we obtain 



4 y/3 5 5 Det'Aj , 
Z[b VDeTA^ 6 

= 0.0047 x (AG) 5 /io ^ exp . (5.43) 

To our knowledge, this formula has been obtained here for the first time, since 
in Refs. [27,13] a closed expression for Z[S 4 ] was not achieved. In particular, 
the isometry factor Q\ was not taken into account and the derivative of the 
C-function was not computed. 
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6 Summary 



Our last step is to use the expressions for Z[S 2 x S* 2 ] and Z[S 4 ] in (5.28) and 
(5.43) and insert these into Eq.(2.15) to find the decay rate 

1 [a $Z{S 2 x S 2 ] 
^V3 Z\S*\ 

= 14.338 y/A (GA)- 2 (/i A) W exp (-^) • (6.1) 

This is the final result of our analysis. This formula gives the rate of semiclas- 
sical decay of de Sitter space due to the spontaneous nucleation of black holes. 
This is the leading mode of decay, since classically de Sitter space is stable [28]. 
The numerical coefficient in the formula originates from the fluctuation deter- 
minants evaluated in the C-function scheme. The factor y/X comes from the 
heat bath temperature coefficient in (2.15) and gives T the correct dimension 
of an inverse time. The coefficient (GA)~ 2 arises due to the combined effect of 
the background isometries. The power of /i A contains the effect of rescalings, 
where we have passed again to the dimensionful renormalization parameter (j, . 
Since quantum gravity is non-renormalizable, ji Q remains undetermined, and 
we have nothing to say about this problem. For numerical estimates it is rea- 
sonable to assume that \i Q ~ G. The last factor in the formula is the classical 
term. The formula is obtained in the one-loop approximation, which is good as 
long as the classical term is large compared to the quantum corrections, that 
is for AG <C 1. Under this condition the nucleation rate is exponentially small. 
Notice that since the overall power of A is positive, the quantum corrections 
provide an additional suppression of the transition rate for small A. 

The formula gives the probability of black hole nucleation per unit proper 
time of a freely falling observer in his Hubble region. The latter is the region 
enclosed inside the observer's cosmological horizon. If a black hole is created, 
then it has the radius 1 / y/X and fills the whole Hubble region. This does 
not mean that the whole space will be eaten by a giant black hole, since de 
Sitter spacetime consists of many Hubble regions, whose number grows as the 
universe expands. Some of these regions will contain a black hole but most 
of them will be empty. The black holes are actually born in pairs, where the 
two members of the pair are created at the opposite sides of the 3-space. The 
interesting conclusion is that for GA <C 1, when inflation is 'slow', the rate 
of black hole nucleation is strongly suppressed, but the created black holes 
are large. This can be understood as a consequence of the fact that the black 
holes are made of the energy contained inside the Hubble region. As the size 
of the latter is large for small A, the created black holes are also large. On 
the other hand, if one is allowed to extrapolate the formula for GA ~ 1, when 
inflation is fast, then the created black holes are small, but they are created 
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in abundance. 



One can see that for late times the number of black holes per unit physical 
volume will be constant. Let us choose for de Sitter spacetime the global 
coordinates associated with the freely falling observers: 

d s 2 = -drf + j cosh 2 | // ) <li% . (G.2) 




Here 77 is the (dimensionful) proper time and dVl\ is the volume element of 
the unit 3-sphere. The volume of the global hypersurface of constant rj 

3/2 ^ 3/2 

is V(rj) = 2tt 2 (f) cosh 3 (yf ^) ~ x (l) exp( v / 3A^). The portion of 
contained inside the future event horizon of any observer has the volume 

Vjj — ^ (x) ^ (^ or ^ a ^ e ^ ms * s ^ ne s P a tial Hubble volume. [This quantity 
slightly depends on the choice of the hypersurface. Even though for any given 
observer one has rj — t, which is the time associated with the observer's 
coordinate system, one has Y* n 7^ E t , unless r/ = t = 0, in which case the 

spatial Hubble volume is Vh = tt 2 (|) ]. As a result, the number of Hubble 
volumes on the hypersurface is Nn(i]) = V(i])/Vn. [One has A^h(O) = 2: the 
de Sitter throat consists of two causally disconnected parts belonging to the 
Hubble regions of two antipodal observers [47].] Multiplying Nn(r)) by T gives 
the black hole nucleation rate per , 

^ = |exp(v^A,)r. (6.3) 



Integrating with respect to rj and dividing by V{rf) yields the average volume 
density of created black holes on S^, 

Pe„ = 4r, (6.4) 



which does not depend on 77. 

The subsequent real time evolution of these black holes is an interesting issue. 
Presumably most of them will immediately evaporate, unless A is very small 
and the black holes are large, in which case however the nucleation rate is 
strongly suppressed. It was argued in [9] that this process could dramatically 
change the global structure of de Sitter space. For more information on this 
issue we refer to [10,9,18] and to the papers cited in Ref.[9]. 

The following steps have been essential in our analysis. We have derived 
Eq.(2.15) for the nucleation rate using the thermal properties of de Sitter 
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space. For this we have approximated the partition function for Euclidean 
quantum gravity with A > by the semiclassical contributions of the S 4 and 
S 2 x S 2 instantons, of which the first yields the free energy F in the Hubble 
volume while the contribution of the second can be regarded as a purely imag- 
inary part of F . In a sense one can think of the created black holes as being 
the bubbles of the new phase spontaneously created out of thermal fluctua- 
tions via quantum tunneling. We have argued that these bubbles may have 
temperature different from that of the heat bath, since they cannot thermalize 
via interactions with the whole reservoir and only exchange energy inside the 
Hubble region. 

To compute the one- loop contributions of the S 4 and S 2 x S 2 instantons we 
have used the standard Faddeev-Popov approach to the path integral. We 
have worked with a one-parameter family of covariant background gauges and 
employed the Hodge decomposition of the fluctuations with their subsequent 
spectral expansion. In our treatment of the conformal modes we have followed 
the standard recipe of complex rotation, up to several lowest lying modes 
for which a different prescription has been applied. In order to integrate over 
zero modes of the Faddeev-Popov operator arising due to the background 
isometries, we have gone beyond the perturbation theory and showed that 
the corresponding integration measure is the Haar measure on the isometry 
group. There are no other zero modes in the problem - for example, the 
standard rotational zero modes are absent because rotations are isometries of 
the backgrounds under consideration. 

We have explicitly determined the spectra of the fluctuation operators. For 
fluctuations around the S 2 x S 2 instanton the spectrum was obtained by 
directly solving the differential equations, while in the S 4 case group theo- 
retic methods have been applied, in which we followed the approach of [27]. 
These spectra have been used in order to compute the functional determinants 
within the ^-function regularization scheme, the corresponding (^-functions be- 
ing studied in detail in the Appendix below. We have checked that our results 
agree with the general formulas for the anomalous scaling behaviour. Finally, 
we have obtained in (5.28), (5.43) the one-loop partition functions for fluctu- 
ations around the S 4 and S 2 x S 2 backgrounds. To our knowledge, in both 
cases such closed expressions have been obtained for the first time. The last 
step has been to use the resulting partition functions in order to calculate the 
nucleation rate T. This describes a constant density of created black holes per 
unit physical volume of the expanding 3-space. 

After the work of Gross, Perry and Jaffe [30] , our analysis presents the second 
example of a complete one-loop computation on a non-trivial background. [f] 

3 Note also that the analysis in [30] was not quite complete, since the spectrum 
is unknown and the (-functions have not been computed, even though the unde- 
termined quantities can be absorbed into the renormalization parameter. We also 
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One may hope that our results can lend further support to the Euclidean 
approach to quantum gravity. 
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| 329 571 

proportional to /x 45 instead of fi 45 [52]. However, since these zero modes are 
in the path integration measure, they do contribute to the anomalous scaling 
on equal footing with all other modes. In fact, the example of flat space gauge 
theories [45] shows that the background symmetry zero modes, when treated 
non-perturbatively as was done above, are of vital importance for obtaining 
the correct running behaviour of the coupling constant. Our result for the 
scaling behaviour agrees with that of Christensen and Duff [13] and with the 
general analysis of Fradkin and Tseytlin [54] . 



Appendix. Calculation of (^-functions. 



In this Appendix we shall study the ^-function 



oo oo 



„ , . ^ v-. (2n + 1) (2m + 1) 

Z(k,L>\s)=>> — ^— 7^ +■ tj, (A.l) 

{(2n+ !)* + (2m + !)* + „}" 1 ' 



which is used in the main text for computing the one-loop fluctuation term on 
the S 2 x S 2 instanton background. Here v is real while A; is a positive integer 

do not understand their treatment of the background isometries and that of the 
non-normalizable deformations of the instanton. 
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such that 2(2k + l) 2 + v > 0. It is assumed that 3?(s) is positive and large 
enough to ensure the convergence of the series. Despite its apparent simplicity, 
the analysis of this ^-function is lacking in the literature. This is probably due 
to the fact that the summation in (A.l) cannot be extended to all integers and 
the standard Poisson resummation techniques do not apply. For this reason 
we use other methods, which are unfortunately rather lengthy. However we 
think that it is necessary to describe the basic steps, especially in view of 
other possible applications of our results. 

In what follows we shall perform the analytic continuation by finding the 
integral representation for Z(k, v\s) that is valid for any s. This will be used 
to compute the values of Z(k,v\0) and j-Z(k,v\s) at s = 0. As a first step, 
we shall consider the related ^-function: 



with (2k + l) 2 + v > 0. The integral representation for this function will be 
useful. In addition, we shall study the (^-function 

h 0X7+3)+"}' 



where k(k + 3) + v > 0, and shall find its value and its s-derivative at s = 0. 
This function is needed in the analysis of fluctuations around the S 4 instanton. 



A.l Computation of Z(k,v\0) and ((k, u\0). 



First we shall compute the values of these functions at s = using the stan- 
dard heat kernel technique. These values determine the scaling properties of 
the system. Later we shall rederive the same values by using the integral rep- 
resentations for Z(k,u\s) and ((k,u\s), and this will provide us with a good 
consistency check. For Q(k,u,c\s) we shall consider only the integral repre- 
sentation, since the values of Q(k,v,c\0) have been computed in [13]. 

A (^-function related to a second order elliptic operator with a positive spec- 
trum can be expressed as 

oo 

as) = —Jt s - l Q(t)dt. (A.4) 
^ ' o 

On compact spaces the heat kernel Q(t) vanishes exponentially fast for large 
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t, while for small t there is the asymptotic expansion 

e(*)~X)c P f, (A.5) 



with r assuming in general both integer and half-integer values. It is not 
difficult to see that 

C(0) = C Q . (A.6) 

The problem therefore reduces to determining the asymptotic expansion of 
the heat kernel. The heat kernels in our problem are given by 

e(k,v\t) = (m-tmy*-'* (a.7) 

for Z(k, u\s) and 

6(k,v\t) = (6(t)-ak\t))e~ ut (A.8) 

for £(k,v\s), where 

oo 

6{t) = £ (2n + l)e-* (2n+1)2 (A.9) 

n=0 

and 

fc-i 

£(k\t) = ^ (2n + 1) e -*( 2n+1 ) 2 . (A.10) 

71=0 

The only difficulty is to find the asymptotic expansions for small t for the 
function 9(t) in (A. 9)0. 6(t) is a partition function for a two-dimensional 
rotator at temperature 1/t. We wish therefore to find its high-temperature 
expansion, and for this we shall construct the integral representation for 9{t). 

Let us consider the "generating function" 

oo 

X ( t ,a) = Y^ e -t^n+lf + i a {2n+l) (An) 
n=0 

4 We note that 9(t) cannot be expressed in terms of theta- functions in a simple 
way, and that the Poisson resummation formula does not directly apply. 
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such that 



9(t) = -ihm^-x{t,a). (A.12) 



x(t, a) fulfills the differential equation 

dx = d^x 
dt da 2 

This has the special solution 



(A.13) 



with the property x(0, a) = 5(a — a ), which allows us to represent the general 
solution of (A.13) as 

oo 

X(t,a)= J x(t,a o )x(0,a o )da o . (A.15) 



The initial value x(0, a ) is obtained directly from the definition (A. 11): 

oo • 

x(0,a o ) = E eia ° (2n+1) = W 2 — . (A.16) 
n=o 2sma 

where we assume that ao has a small positive imaginary part in order to ensure 
convergence of the geometrical series. We can now insert this into (A.15) and 
the result into (A.12). Introducing the new variable x = a^/A we obtain the 
sought for integral representation 

6(t) = -L= Je-*'* dX ^ . (A.17) 
W v& J sin(2v^) 1 ; 



Here we should remember that x has a small imaginary part, such that the 
integration is actually performed along a contour parallel to the positive real 
axis and approaching it from above. 

It is now a straightforward task to find the asymptotic expansion of the integral 
in (A.17) for small t, since the only non-trivial contribution comes from a small 
neighbourhood of x = 0: 

0(t)~l (l + it + ^t 2 + 0(t 3 )). (A.18) 
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Inserting this into (A. 7) and (A. 8) gives the asymptotic expansions for the heat 
kernels Q(k, v\t) and 9(k, v\t), whose coefficients C determine the ^-functions 
at s = 0: 



Z(k,u\0) 




360 ' 



13 



(A.19) 



and 



ak,v\o) = ±-±v-k 2 . 



(A.20) 



To check these results we note that the definitions in (A.l) and (A. 2) imply 
that 



with fc 2 > k\. Setting here s = we obtain a non-trivial relation for Z(k, u\0) 
and £(k,is\0), and this is fulfilled by the expressions in (A.19) and (A.20). 
Finally we use (A.19), (A.20) to obtain the values used in the main text: 



Z(2,-10|0) = ^-, C(2, -9|0) 1 , C(2, -1|0) = -— . (A.22) 



A. 2 Computation of ((k,u\s), and Q(k, v, c\s). 

It is usually more difficult to determine the derivative of a ^-function at s = 
than the value of the function itself, since the knowledge of its behaviour in a 
neighbourhood of s = is required. We shall perform the analytic continuation 
of the C-functions defined by Eqs. (A.l), (A. 2) and (A. 3) to arbitrary values 
of s with the use of the relation sometimes called 

A. 2.1 The Abel- Plan formula. 

This can be derived using the obvious relation 



k 2 -l 



Z(h, u\s) = Z(k 2 , v\s) + 2 Y, ( 2m + !)C(fe, v + (2m + l) 2 |s) 




{(2n+l) 2 + (2m+l) 2 + z/} 



(2n + 1) (2m + 1) 



(A.21) 




(A.23) 
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where the contour C encompasses the part of the real axis with Re(z) > k 
(see Fig. 4) and f(z) is analytic for Re(z) > k. The idea is to split C into three 
parts, C\ + C 2 + C 3 , as shown in Fig. 4. For the first part, C±, the integral can 
be written as 

/ (id^ - / W = //(*) * + / ^ > (A- 24 ) 



where in the integral over C\ on the right the contour is then rotated to the 
position C\ as shown in Fig. 4. Such a rotation is possible if only f(z) tends 
to zero fast enough for Re(z) > k and \z\ — > oo. 




-X X X X X- 



X X X X X X- 



Z+(T„) 



Fig. 4. Left: Starting from the contour C\ + C2 + C3 and rotating we arrive 
at C\ + C2 + C3. Right: the same when a branching point at z = z + {t) is 
present. The contour C\ will then wrap around the cut leading to the additional 
contribution due to C. The point z = z + (t) is in the region of interest for 

T > T*. 

The integral over the second portion of the contour, C2, is equal to ^f(k), 
while in the integral over C 3 the contour is rotated to the position C 3 as 
shown in Fig. 4. As a result, we arrive at the Abel-Plan formula 



n=k 



- 2 f(k) 



f(t)dt + i J 




f(k + it)-f(k-it) 



dt. (A.25) 



This formula can be used for analytic continuation of C-fimctions, in which 
case f(t) depends also on s, f — f(t,s). The analytic continuation to small 
values of s is performed in the first integral on the right in (A.25). This usually 
converges only for 9?(s) large and positive, but can often be computed in a 
closed form, and then one can continue the result to arbitrary s. The second 
integral on the right in (A.25) usually cannot be computed in a closed form, 
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but it converges for any s. Let us first apply the Abel-Plan formula to the 
^-functions in (A. 2) and (A. 3). 

A. 2. 2 Analytic continuation of ((k,u\s). 

Applying (A. 25) to the series for ((k, z/|s) in (A. 2), we have 

which is analytic for ^St(z) > —1/2 and decays fast enough for \z\ — > oo 
provided that 9ft(s) is large enough. As a result, we can use the Abel-Plan 
formula, which gives 



((k,v\s) = [k 



2 J {(2k + l) 2 + u} s 4(s — 1) {(2k + l) 2 + i/} s_1 

idt ( 2k + 1 + lit 2k + l-2it \ 

e 2 -* - 1 \{(2k + 1 + 2it) 2 + u} s ~ {(2k + 1 - 2it) 2 + u} s ) ' ( " 



This representation is finite for all s, apart from s — 1, where the pole is 
located. The remaining integral here converges uniformly for \s\ < oo, which 
allows us to differentiate with respect to s. If we set s = 0, then the integral 
can be easily computed. We find ((k,u\0) = j^ — ju — k 2 , and this agrees with 
the value obtained above in (A. 20). 

Similarly, we can differentiate (A. 27) with respect to s and then set s = 0. 
This gives 



('(k, u\0) = ^ W(\n W - 1) - (k + In W 

oo 

r dt 

+2j ^ f — [ (t\nA+(2k+m, (A.28) 
o 

where W = (2k + l) 2 + v and 



A=(W- At 2 ) 2 + 16(2A; + l) 2 t 2 , * = arctan A{2k + l \ t . ( A .29) 

W - At 2 



For any k and v the integral in (A.28) is convergent and can be evaluated 
numerically. Notice that ('(k, v\0) is not needed in the main body of the paper, 
and for this reason we do not quote the actual number here. 



55 



A. 2. 3 Analytic continuation of Q{k,u, c\s) . 

The procedure is exactly the same as above. Denoting 



= (2, + 3)(,(» + 3) +C ) 

{z(z + 3) + v} s v ' 



the direct application of the Abel-Plan formula (A. 25) yields 



W 2 ~ s 
2 



Q(k, v, c\s) = (k + ^ (k(k + 3) + c) W~ s + -i- 

oo 

+^ W 1 "' + / <(/(* + «) " /(* - «)) , (A-31) 



where W = k(k + 3) + v. Setting s = the integral can be easily computed 
leading to 

Q(k, v, c|0) = A; 4 - 2A; 3 - (c + ^) ^ 2 (A.32) 
+ (3-2cH + ^ + g-,)c-±±. 
Next, differentiating (A. 30) with respect to s and setting s = gives 



zy, c|0) = + (W + c-*/) mW 
+ - (in W - ^ W 2 + (c - i/) (In W - f ) W + Q . (A.33) 



Here 



oo 

dt 



/at 
27rt _ {t (6fc(fc + 3) + 2c + 9 - 2t 2 ) In .A 

o 

+ (4A; 3 + 18A; 2 + (f 8 + 4c) k + 6c - 6 (2ib + 3) t 2 )^} (A.34) 

with 

A = t i + (2W-4u + 9)t 2 + W 2 , * = arctan ( 2fc + 3 )* . (A.35) 

W — t 2 

Evaluating the integral numerically, the two values used in the main text are 
Q'(2,0,-4|0) = 3.72344, Q'(2, -4, 0|0) = 6.65246 . (A.36) 
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Finally, for the function ((s) = 3 S (^Q(2, 0, — 4|s) - \Q(2, -4, 0|s)) used in 
Eq.(5.41) in the main text one obtains with the help of (A. 32) and (A. 36) 

C(0) = w , C'(0) = T 1 = 6.10158. (A.37) 
A. 3 Computation of Z(k,u\s) 

Let us not turn to our main task - the evaluation of the double-sum function 
Z(k, v\s), which has been defined for large values of 9R(s) by (A.l). The idea 
is to express it in terms of the single-sum function ((k,u\s). 

It follows from the definitions (A.l) and (A. 2) that 

oo 

Z(k, is\s) = J2 ( 2n + l )C(k, v+(2n + lf\s) . (A.38) 

n=k 

Here we can use the integral representation (A. 27) for ((k, v + (2n + l) 2 |s). 
Indeed, if v is real and (2k + l) 2 + v > then the same remains true upon 
replacement v — > u+ (2n+ 1) 2 , and the formula (A. 27) therefore applies. Now, 
replacing in (A. 27) v by v + (2n + l) 2 and assuming for a moment that 3?(s) 
is large and positive, the integral in (A. 27) converges uniformly with respect 
to n for n — > oo. This allows us, upon insertion of (A. 27) into (A.38), to 
interchange the orders of summation and integration. The result then can be 
extended to any s by analytic continuation. This gives 



/ 1 \ 



2n+ 1 

„ k {(2k + iy + v + (2n + iyy 



1 ~ 2n + 1 



4(s - 1) £ i( 2k + !) 2 + » + (2" + I) 2 } 5 " 1 

2n + 1 

{(2A; + 1 + 2ir) 2 + z/ + (2n + 1) 2 } S 



1 idr ( 

n \ n=/e 



- (a+i - aT) | { (2. +1 -2^:: +(2n+1 p } i < a -' 

One can see that all sums here are exactly the same as in the definition of 
C(/c, f|s) in (A. 2) - up to the replacements v — > z/ + (2fc+ l) 2 and z/ — > z/(r) = 
z/ + (2/c + 1 + 2ir) 2 . Since the definition in (A. 2) makes sense for arbitrary 
values of v (the series always converges for 9?(s) big enough), we can express 
the sums in (A. 39) in terms of ((k, u+ (2k + l) 2 |s) and ((k, u(t)\s). This leads 
to the the following formula: 
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Z(k, u\s) = { k + \) C(k, v+(2k + l) 2 \s) (A.40) 
-C(k, v + (2k + lf\s - 1) + J -^— x {T(r) ~ H-r)} , 



4(s 

with T(t) = (2k + 1 + 2%r) ((k, u(t)\s). In this formula the first two terms on 
the right are determined by the integral representation (A. 27) for arbitrary s. 
We are left with computing the remaining integral over r. The problem here 
is that the parameter v(r) is complex, and for this reason we cannot directly 
apply the integral representation (A. 27) to compute ((k, v(t)\s). 

Let us recall that the formula (A. 27) was derived assuming that the function 
f(z) in Eq.(A.26) had no poles for ^St(z) > k. This allowed us to rotate the 
integration contour as shown in the left part of Fig.4 without intersecting 
singularities. Let us now replace v by v(r) = v + (2k + 1 + 2ir) 2 . As a result, 
f(z) in Eq.(A.26) is replaced by 

v _ 2,2 + 1 2z + 1 

m ~ {(2z+1)* + u(t)Y ~ {4(z - z + (t))(z - z.(t))}" 1 ' > 



with z±(t) = \(-l ± i^Ar)). For r = one has 3£(z±(0)) = -\. As r 
increases, the point z + (r) moves to the right in the complex plane (while Z-(t) 
moves to the left), but as long as $1(z + (t)) < k one can still use the formula 
(A. 27). However, for large enough values of r the pole at z — z + (r) enters the 
region of interest, that is the part of the complex plane with $l(z) > k, and 
we can no longer use the formula (A. 27). 

To tackle the problem we notice that the pole of f(z) at z = z+(t) is a 
branching point, and one can choose the cut in the complex plane as shown 
in the right part of Fig.4. We then repeat the steps leading to the Abel- 
Plan formula and the additional problem we encounter is the following: when 
we rotate the integration contour as we did before, the contour will wrap 
around the cut as shown in Fig.4. The resulting contour will then consist of 
two disconnected pieces. The first piece will be the same as the old contour 
C\ + C2 + C3 (see Fig.4). The second piece is the contour C wrapping around 
the cut. Integrating around such a combined contour, the result will consist 
of two parts, 

((k,v(T)\s) = ( old (k,v(T)\s) + 9(T-n) J J^dz. (A.42) 

c 



Here the first term on the right, Coid(&> h/ { T )\ s )i is the function given by the 
previous expression in (A. 27) with v being replaced by v(j). The second term, 
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with f(z) given by (A. 41) and the contour C as shown in the right part of 
Fig. 4, is the contribution of the cut. The step function 9(t — r*) reflects the 
fact that the cut contributes only for large enough r when the pole enters 
the region $l(z) > k. Here 6{x) = for x < and 9{x) = 1 for x > 0, and 
U(z+(t*)) = k. 

The representation (A. 42) applies for all values of s and for any r > 0. Sim- 
ilarly, one can obtain ^(k,v(—r)\s) (the cut then resides in the upper half- 
plane). As a result, the function F{t) — F(—t) in the integrand in Eq.(A.40) 
is defined for any r > 0, and the integral converges due to the damping expo- 
nential factor. This finally gives Z(k, u\s) for any s. 

Let us first check our result by computing Z(k,u\0). For s = the function 
f(z) has no poles and the contribution of the cut vanishes. The remaining 
integrals then can be easily computed, which gives for Z(k,v\0) exactly the 
same expression as in Eq.(A.19). 

Let us now compute Z'(k,u\0). Since all integrals in (A. 40), (A. 41) converge 
uniformly with respect to s (at least for \s\ < oo), we can differentiate the 
integrands with respect to s and then set s = 0. The result can be represented 
in the following form: 



oo 

z'(kMo) = n + 2f^±-g(t) 



oo oo 





Here 



H = (-2k A + (2 - -) k 2 + k + — (4 + 4i/ - v 2 )) ln( 2 (2k + l) 2 + u) 
V 2 32 / 

+3A; 4 + 2A; 3 + ^ [v - 2) k 2 + - [v - 6) k + ^ (3z/ 2 - 4i/ - 20) . (A.44) 



In addition, 



Q(t) = l - (At 2 - lQk 2 - 12k - 2 - v) In P (A.45) 
+ (2k + l)(6t 2 - 2k (2k + 1) - V -) $ , 

where we have used 

P = u 2 + (4(2Jfe + l) 2 - 8t 2 ) v + 4(2A; + l) 4 + 16t 4 , 
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$ = Phase[(2A; + l) 2 + v/2 - 2t 2 + i 2(2k + l)t\ , (A.46) 
and —n <Phase[x + iy] < ir is the phase of the complex number. Next, 

W(r,t) = {((2fc + l) 2 -AtT)\nQ (A.47) 
+ 4 (t - r)(2k + 1)*} - {(t, t) <-> (t, -r)} (A.48) 

with 



g = 16 (t 2 + r 2 ) 2 + (v + 2) 2 - 81/ [t 2 + r 2 ) + 128 {k 2 + k + -) t r 

+ 16 fc(ife + 1) i/ + 32 k(k + l)(2k 2 + 2k+l), 
* = Phase[(2A; + l) 2 + v/2-2 (t 2 + r 2 )+i2 (2k + l)(t- r)] . (A.49) 



Finally, the contribution of the cut is 



oo , oo 



r* 



where z(r) = -| + ^r 2 - (2ib + l) 2 - v + 4i(2A; + l)r, and 3?(z(r*)) = fc. 

We now use the formulas above in order to evaluate Z'(2, — 10 1 0) , which value 
is needed in the main text. Setting k = 2 and v = —10 we obtain for the 
first term on the right in (A. 43) H. = 1.9445. The second term, containing 
the integral over t, is evaluated numerically to give —19.9469. The numerical 
value of the term containing the double integral is —0.1294. As for the last 
term, S, it is exponentially small and is of the order of 10~ 12 . This is because, 
as one can see from (A.50), the value of S is suppressed by the factor of 
exp{-2n(T* + S(z(t*))} = exp{-47rv / 5}. 

Summing everything up, we obtain 

T = Z'(2, — 10|0) = -18.3118 . (A.51) 
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